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Chapter  1  Introduction 


Section  1-1  Introduction 

During  the  1980’s,  advances  in  computer  hardware  and  software  have  made 
it  possible  to  perform  ab  initio  calculations  on  medium  sized  molecules  at  levels  of 
sophistication  that  are  capable  of  yielding  chemically  useful  and  reliable  results.1 
Indeed,  ab  initio  calculations  offer  a  probe  into  molecular  behavior  that  may  not  be 
easily  accessible  with  experimental  techniques.  Ab  initio  calculations  may  be 
regarded  as  another  tool,  equal  to  experiment,  to  gain  further  insight  into  chemical 
systems. 

Using  ab  initio  calculations  to  supplement  existing  experimental  data,  I  have 
investigated  a  variety  of  chemical  phenomena  These  include:  the  effects  of 
substituents  on  altering  carbon-hydrogen  and  silicon-hydrogen  bond  strengths,  the 
mode  by  which  the  empty  orbital  on  boron  alters  the  physical  properties  of  an 
adjacent  phosphorous  atom,  and  the  effect  of  substituents  on  the  singlet-triplet  energy 
gap  in  a  non-Kekul6  molecule. 

Each  of  these  topics  is  covered  separately  in  the  following  chapters;  but  first, 
a  brief  discussion  of  ab  initio  theory  is  presented,  in  order  to  familiarize  the  reader 
with  the  computational  methodology  that  has  been  used. 
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Section  1-2  Overview  of  Ab  Initio  Theory 

Quantum  chemistry  is  comparatively  young,  having  its  origins  in  the  mid-1920’s 
with  the  introduction  of  Schrodinger’s  equation.  This  equation  describes  the 
behavior  of  particles,  such  as  electrons  and  nuclei,  which  comprise  atoms  and 
molecules.  Thus,  at  least  in  principle,  the  physical  solutions  of  the  Schrodinger 
equation  should  allow  the  prediction  of  physical  properties  and  chemical  reactions 
of  atoms  and  molecules.  First,  however,  methods  had  to  be  found  to  obtain  a 
mathematical  solution  to  the  time  independent  Schrodinger  equation. 


Section  1-2-1  The  Time  Independent  Schrodinger  Equation 
The  time  independent  Schrodinger  equation  is: 

HH  -  EH 

where  A  is  the  Hamiltonian,  an  operator,  which  defines  the  allowed  wave  functions, 
Y,  and  associated  energy  levels,  E,  for  a  particular  system.  The  Hamiltonian  is 
comprised  of  two  distinct  terms,  one  describing  the  potential  energy  of  the  system, 
V,  and  another  for  the  kinetic  energy  of  the  system,  * . 

H  -  V  +  f  (2) 


all  particles 

>>-£  E 

i  i<J 


flL L) 

rv  ) 


(3) 


3 


.  2  all  particles  , 

rp  n  1 

‘  ’  "2  r  ^ 


2  \ 


jl  +  _ai  +  _a_ 

axf  +  by)  +  az? 


(4) 


Where  e(  is  the  charge  of  each  particle,  nij  is  its  mass,  and  is  Planck’s  constant 
divided  by  2 it. 

Since  the  mass  of  an  electron  is  about  one  thousandth  that  of  a  proton,  the 
electrons  should  be  able  to  adjust  their  motion  rapidly  to  adjust  to  any  changes  in 
position  of  the  nuclei.  Hence,  separating  the  motion  of  the  electrons  from  that  of 
the  nuclei  is  possible  and  makes  solving  the  Schrodinger  equation  much  easier.  This 
separation  of  motion  is  referred  to  as  the  Born-Oppenheimer  approximation.2 

This  approximation  allows  transformation  of  the  Schrodinger  equation  to: 


HeUc(r,R)  Y<fec(r)  -  Etkc(R)  TeUc(r) 


(5) 


where  the  electronic  Hamiltonian  depends  on  both  the  positions  of  the  electrons  (r) 
and  the  nuclei  (R),  but  operates  only  on  the  wave  function  for  the  electrons.  For  a 
particular  set  of  R,  the  electronic  wave  function  is  an  explicit  function  only  of  the 
positions,  r,  of  the  electrons.  Solutions  of  equation  (5)  thus  give  the  electronic  wave 
functions  and  associated  energies  at  a  given  set  of  values  of  R.  Nuclear  motion  can 
be  determined  by  solutions  of  a  nuclear  Schrodinger  equation  in  which  the  nuclei 
move  in  an  effective  potential,  which  is  provided  by  the  electrons  and  is  equal  to 
P'*C(R). 
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Section  1-2-2  Potential  Energy  Surfaces 

The  electronic  wave  function  of  most  interest  is  usually  the  lowest  energy 
solution  of  equation  (5),  and  this  wave  function  is  called  the  ground  state  wave 
function. 

The  variations  with  R  of  the  energy,  Flec(R),  including  the  nuclear  repulsion 
energy,  for  the  ground  state  wave  function  gives  the  ground  state  potential  energy 
surface.  Stationary  points  on  this  surface  represent  points  that  have  zero  first 
derivatives  of  the  energy  with  respect  to  all  nuclear  motions.  Hence,  at  stationary 
points  the  nuclei  have  no  net  force  exerting  on  them. 

If  a  stationary  point  is  an  energy  minimum,  all  the  second  derivatives  of  the 
energy  at  this  point  are  positive.  Such  a  point  may  only  be  a  local  minimum.3  The 
global  minimum  --  the  minimum  with  the  lowest  energy  on  a  surface  -  is  commonly 
referred  to  by  chemists  as  the  equilibrium  geometry. 

Another  type  of  stationary  point  of  interest  on  the  potential  surface  is  a 
transition  state.  Transition  states  are  saddle  points,  i.e.  they  are  stationary  points  at 
which  the  energy  is  a  maximum  for  one  direction  of  nuclear  motion.4 
Mathematically,  the  matrix  of  second  derivatives  of  the  potential  surface  -  the 
Hessian  matrix  --  has  exactly  one  negative  eigenvalue.  Transition  states  are  the 
highest  energy  point  along  the  lowest  energy  pathway  connecting  two  minima. 
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Section  1-2-3  Expressing  the  Wave  Function  as  Molecular  Orbitals 

For  one-electron  atoms  or  molecules,  exact  solutions  to  equation  (5)  can  be 
obtained  in  closed  form.  For  many-electron  (>  1)  systems,  this  is  not  possible  and 
the  first  step  in  solving  equation  (5)  is  determining  an  appropriate  representation  of 
the  electronic  wave  function.  If  electrons  are  assumed  to  behave  as  independent 
particles,  each  electron’s  spatial  wave  function  motion  is  described  with  its  own 
orbital,  f.  The  square  modulus,  |  4r  |2,  represents  the  probability  distribution  of  the 
position  of  an  electron  in  tjr  spatially. 

In  addition  to  a  spatial  wave  function,  each  electron  also  has  a  spin  wave 
function  with  a  spin  quantum  number  of  ±  %.  The  two  possible  spin  wave  functions 
are  usually  referred  to  as  a,  and  /?.  As  a  result,  the  space-spin  orbital  for  a  single 
electron  is  either  <Jra(x,y,z)  or  \|i*(x,y,z).  These  space-spin  wave  functions  are 
sometimes  abbreviated  as  x„(i)  where  i  identifies  which  electron  occupies  this 
space-spin  orbital,  Zn- 

If  electrons  are  approximated  to  behave  as  independent  particles,  one  might 
be  tempted  to  write  the  overall  wave  function,  7,  as  a  simple  product  of  the 
individual  space-spin  orbitals: 

7  -  X,(l)  Xj(2)  -  X„(")  (<■) 

However,  since  electrons  are  indistinguishable,  fermion  particles,  if  any  two  electrons 
are  transposed,  the  sign  of  the  wave  function  must  change.  Hence,  the  wave  function 
must  be  antisymmetrized.  This  is  most  easily  accomplished  by  writing  the  wave 
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X,(D  x2(D  •••  X,(D 

^  Xi(2)  X2(2)  •  X„(2)  (7) 

Xi(«)  X2(«)  X„(«) 

function  as  a  Slater  determinant,  as  expressed  in  equation  (7).  Since  a  determinant 
goes  to  zero  if  any  two  columns  are  the  same,  it  is  clear  that  each  electron  must  have 
a  unique  %n •  This  gives  rise  to  the  Pauli  exclusion  principle6  which  states  that  two 
electrons  cannot  have  the  same  spin  and  space  wave  functions. 

Writing  Y  as  a  determinant  is  rather  cumbersome.  Y  can  be  expressed  more 
conveniently  in  Dirac  notation  as  either  a  bra  or  a  ket,  using  the  space-spin  functions 
found  on  the  principle  diagonal  of  the  Slater  determinant.  Thus,  in  this  shorthand 
designation,  Y,  written  as  a  ket,  can  be  abbreviated 

T  -  I  X,(D  X2m  X»  )  (8) 

Next  to  be  determined  is  an  appropriate  mathematical  expression  for  the 
spatial  part,  i|r,  of  each  %  in  Y.  If  Y  were  the  wave  function  for  a  single  atom,  atomic 
orbitals,  <p,  such  as  Slater-type  atomic  orbitals,6  would  be  appropriate  for  Y-  For 
molecules  comprised  of  many  atoms,  an  approximate  method  is  to  describe  each 
molecular  orbital  (MO),  ijrn,  as  a  linear  combination  of  atomic  orbitals  (LCAO).7 
Mathematically, 

♦  »  -  E 

i 


(9) 
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where  each  i|rn  has  its  own  unique  and  orthogonal  set  of  cni,  but  all  i|rn  are  formed 
from  the  same  atomic  orbitals  <px.  Hence,  for  a  given  geometry,  any  two  unique  MOs 
will  differ  only  by  the  coefficients  representing  the  contributions  of  each  atomic 
orbital. 

Section  1-2-4  Hartree-Fock  Theory 

The  electronic  wave  function  of  interest  most  often  is  the  ground  state  wave 
function.  One  way  to  find  the  lowest  energy  wave  function  is  by  application  of  the 
variation  principle.8  Assuming  that  Y  has  the  form  in  (7),  the  spatial  parts,  of  Y 
are  varied,  subject  to  the  constraint  that  the  i|r  remain  orthogonal  and  normalized. 

0°) 

This  gives  rise  to  the  Hartree-Fock  equation^’9  for  finding  the  optimal  to  use  in 
the  wave  function,  Y,  for  the  ground  state, 

Hi  ~  e,t, 

where  F  is  called  the  Fock  operator  and  e{  is  the  energy  of  orbital,  ty. 

Although  the  Fock  operator  is  a  one-electron  operator,  it  includes  the  effect 
of  Coulomb  repulsion  between  electrons  by  containing  terms  for  the  effective  field 
felt  by  each  electron.  Thus,  in  order  to  compute  f%  one  must  already  know  Y.  In 
practice,  this  involves  using  successive  approximations  -  guessing  Y,  using  this  guess 


8 

to  compute  I7,  finding  a  new  set  of  orbitals,  i|r,  with  which  to  construct  ¥,  recomputing 
F,  and  continuing  on  until  the  process  converges  to  a  self-consistent  solution.  Thus, 
Hartree-Fock  (HF)  theory  is  sometimes  also  called  Self-Consistent  Field  (SCF) 
theory. 

If  the  LCAO  approximation  in  equation  (9)  is  made,  the  energy,  E,  must  be 
minimized  with  respect  to  variations  in  the  coefficients  cni.  The  coefficients  of  the 
atomic  orbitals  in  the  LCAO-MOs  must  again  satisfy  the  orthonormalization 
conditions, 

/*!*/*  “  EE  CniSHmCmj  -  6 ij  (12) 

*'l  *1-1 


where  S^m  is  the  overlap  matrix  between  atomic  orbitals,  whose  elements  are  defined 
by 

-  jV-m-d)*  <13> 


The  LCAO  approximation  leads  to  the  Roothaan-Hall  equation 


W-l 


-0 


(14) 


for  finding  the  coefficients  in  the  MO’s  ty,  that  minimize  the  total  energy.  In  the 
above  equation,  €;  is  again  the  one-electron  energy  of  LCAO-MO  ty;  and  Fnm  is 
called  the  Fock  matrix,  whose  elements  are  defined  by 
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where  FFore(l)  is  comprised  of  those  terms  in  the  Hamiltonian  that  operate  on  one 
electron  at  time,  kinetic  energy  and  nuclear-electron  attraction. 


(1) - III  +  S—  + 


Thus,  the  first  term  in  equation  (15)  gives  the  kinetic  energy  of  and  the  nuclear 
attraction  energy  experienced  by  the  electron  in  the  region  where  <pn  and  <pm  overlap. 
The  second  term  in  equation  (15)  represents  the  average  Coulombic  repulsion 
experienced  by  an  electron  in  this  region,  where 

(nm\Xa)-fJ(t>l( l)<l>m(l)  — W(2)4>0(2)  dxldytdzldx2dy2dz2  (18) 

\ri2) 

is  the  classic  Coulomb  repulsion  energy,  and  the  third  term, 


(nXlmo)-J|<j>*B(l)<t>1(l)  —  4>*w(2)<{>0(2)  dx^dzfady^ 

Viz) 


(19) 
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is  a  correction  that  comes  from  antisymmetrization  of  the  electronic  wave  function 
(Pauli  exclusion  principle)  and  from  the  fact  that  electrons  do  not  repel  themselves. 
Finally, 


occ 


*  2  E  C^mi 

i-1 


(20) 


is  the  bond  order  between  atomic  orbitals  <pn  and  <£m. 

The  total  energy,  E,  of  a  system  can  be  calculated,  using  components  of  the 
Roothaan-Hal!  equations,  by  adding  the  nuclear  repulsion  energy  to  the  purely 
electronic  energy. 


N  N 


M  A- 1  7  7 


i»-l  m-1  /<-l  U-l  KAB 


(21) 


As  with  the  Hartree-Fock  equation,  the  main  difficulty  in  solving  the 
Roothaan-Hall  equation  (14)  is  that  the  Fock  matrix  is  dependent  on  Y,  through  the 
dependence  of  Fnm  on  the  bond  order,  Pnm,  in  equation  (15).  Therefore,  solving 
equation  (14)  is  an  iterative  process.  A  guess  is  made  for  Y,  which  is  then  used  to 
calculate  Pnm,  which  is  required  to  compute  each  element,  Fnm,  of  the  Fock  matrix. 
Then  solving  equation  (14)  for  cni  gives  a  better  guess  for  Pnm  which  is  then  used  in 
turn  to  determine  a  new  Fock  matrix.  This  cycle  is  repeated  until  self-consistency  is 


reached. 
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Section  1-2-5  Closed  Shell  Systems 

In  general,  closed  shell  systems  are  easier  to  calculate  than  open  shell  systems. 
With  a  closed  shell  system,  two  electrons  occupy  each  MO;  so  there  is  no  need  to 
differentiate  the  orbitals  occupied  by  the  a  spin  electrons  from  those  occupied  by  the 
0  spin  electrons.  The  two  sets  of  MOs  can  have  the  same  spatial  functional  form; 
hence,  only  one  set  of  coefficients  for  the  atomic  orbitals  is  needed  in  order  to 
completely  specify  the  wave  function  for  each  pair  of  electrons.  This  method  is 
designated  Restricted  Hartree-Fock  (RHF)  since  the  orbitals  for  electrons  of 
opposite  spin  are  restricted  to  be  the  same. 

Section  1-2-6  Open  Shell  Systems 

Some  features  of  open-shell  molecules  cannot  be  adequately  described  using 
a  RHF  wave  function.10  According  to  the  Pauli  exclusion  principle,  two  electrons 
of  the  same  spin  cannot  be  in  the  same  region  of  space  at  the  same  time;  hence, 
electrons  of  the  same  spin  are  always  correlated  in  an  antisymmetrized  wave 
function.  However,  the  same  is  not  true  of  electrons  of  opposite  spin.  Since  an 
open-shell  system  has  one  or  more  singly  occupied  orbitals,  the  forces  acting  on  each 
a  spin  electron  are  usually  different  from  those  acting  on  each  /3  spin  electron.  As 
a  result,  the  optimal  MOs  containing  the  a  spin  electrons  are  not  necessarily  spatially 
identical  to  those  containing  the  /3  spin  electrons.  If  the  molecular  orbitals  for 
electrons  of  opposite  spin  are  not  restricted  to  be  identical,  a  self-consistent  field 
calculation  is  referred  to  as  unrestricted  Hartree-Fock  (UHF). 
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Section  1  ^  Beyond  a  Single-Determinant  Wave  Function  Solution 

Approximate  solution,  via  Hartree-Fock  theory,  of  the  time  independent 
Schrodinger  equation  using  a  single-determinant  wave  function  is  sufficient  to 
describe  many  molecular  systems.1  However,  Hartree-Fock  theory  does  not  provide 
a  correct  treatment  of  electron  correlation.  Electrons  of  the  same  spin  are  partially 
correlated  via  antisymmetrization;  but  electrons  of  opposite  spin  are  not  correlated 
at  all  with  a  RHF  wave  function.  When  electron  correlation  is  important,  a 
single-determinant  wave  function  does  not  provide  even  a  qualitatively  adequate 
description  of  the  wave  function;  and  multiple-determinantal  wave  functions  are 
required.1 

Correlation  of  the  electrons  serves  to  reduce  the  probability  of  two  electrons 
simultaneously  occupying  the  same  region;  and,  thus,  correlation  reduces  the 
Coulomb  repulsion  energy.  The  need  for  including  electron  correlation  becomes 
apparent  from  consideration  of  the  dissociation  of  Hj  to  form  two  neutral  hydrogen 
atoms.  The  Hartree-Fock  solution  of  the  wave  function  at  any  separation  is  given 
by 

Y  -  |o“(l)op(2)>  (22> 


where,  in  LCAO-MO  approximation 
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a 


\/2  +  2S 


{*A  +  <M 


(23) 


The  spatial  portion  of  this  wave  function  can  be  expanded  to  give 

*  -  ^[+<(1)*^2)  *  */l)*^(2)  +  4>,(1)*,,(2)  *  <t>^l)<t>»(2)]  (24) 


This  predicts  that  at  any  internuclear  separation  of  atoms  A  and  B,  the  electrons  will 
be  localized  to  separate  hydrogen  atoms  only  half  the  time;  the  other  half,  they  will 
both  be  on  the  same  hydrogen.  However,  at  infinite  separation  of  A  and  B,  the 
electrons  should  always  be  on  different  nuclei.  This  requires  that  the  wave  function 
at  infinite  separation  be 

Y.  -  <t>,(2)  +  4>*(1)  <M2>]  (25) 

v2 

which  cannot  be  written  as  a  single  determinant  using  the  molecular  orbitals  from  the 
Hartree-Fock  wave  function. 

In  order  to  obtain  (25),  a  second  determinant  containing  the  a  anti-bonding 
orbitals  must  be  included.  A  wave  function,  general  enough  to  give  equation  (25)  at 
large  internuclear  distances,  can  be  written,  using  the  Hartree-Fock  orbitals  as 
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Y  -  Cj  |  o*(l)  op(2)  >  -  c2 1  o*  “(1)  o*  p(2)  )  (26) 

where 

(27) 

since  7  of  equation  (26)  can  be  expanded  to 

Y-cc[<|>a(1)<J>,(2)  +  4>*(1)4>X(2)]  +  c,[<|>A(l)<|>A(2)  +  4>B(1)<|>B(2)]  (2g) 

where  the  coefficient,  cc  «  q  +  Cj,  of  the  covalent  terms  and  that,  Cj  «  q  -  Cj,  of  the 
ionic  terms  are  no  longer  restricted  to  contribute  equally. 

The  most  straightforward  way  of  correlating  the  motions  of  many-electrons  is 
to  obtain  variationally  a  multiple-determinant  wave  function,  formed  as  a  linear 
combination  of  Slater  determinants,  which  are  each  comprised  of  orbitals  obtained 
form  Hartree-Fock  theory.  This  treatment  is  called  configuration  interaction  (Cl). 
Alternately,  instead  of  finding  the  wave  function  and  its  energy  variationally, 
perturbation  theory  can  be  applied.  A  third  approach  is  to  start  with  a 
multiple-determinant  wave  function  and  simultaneously  optimize  both  the  mixing 
coefficients  of  the  determinants  and  the  orbitals  in  them.  This  approach  is  called 
Multi-Configuration  Self-Consistent  Field  (MCSCF). 
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The  main  premise  behind  most  Cl  and  perturbation  calculations  is  that  the 
Hartree-Fock  wave  function  is  a  reasonable  approximation  to  the  actual  wave 
function.  Therefore,  it  is  reasonable  to  write 

(29) 

E(exact)  =  E(Hartree-Fock)  +  E(correlation) 

where  E(exact)  is  the  lowest  possible  energy  one  can  obtain  from  a  given  set  of 
atomic  orbitals  after  completely  accounting  for  electron  correlation.  Ideally,  this 
energy  contribution,  E(correlation),  is  small  and,  hence,  as  discussed  later,  may  even 
be  calculated  from  perturbation  theory. 

The  electronic  configurations  whose  mixing  with  the  ground  state  wave 
function  will  lower  the  overall  energy  are  those  of  the  same  spin  multiplicity  and 
overall  symmetry  as  this  wave  function.  These  configurations  are  generated  by 
exciting  one  or  more  electrons  from  the  ground  state  HF  orbitals  to  higher, 
unoccupied,  virtual  orbitals.  Starting  from  HF  orbitals,  Brillouin’s  theorem 
guarantees  that  only  doubly  excited  configurations  mix  directly  with  the  HF 
configuration.11  The  coefficients  of  each  excited  configuration  can  either  be 
calculated  variationally  or  obtained  from  perturbation  theory. 

Section  1-3-1  The  Different  Levels  of  Cl  Calculations  Possible 

In  a  Cl  calculation,  one  must  determine  the  interaction  energy  of  each 
configuration  with  all  the  others,  store  the  resulting  matrix,  and  obtain  the  lowest 
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eigenvalue  and  eigenvector.  Since,  the  number  of  configurations  grows  very  rapidly 
with  the  number  of  electrons  and  orbitals,  one  does  not  always  have  enough  disk 
space  to  store  all  the  Cl  matrix  elements  or  have  enough  CPU  time  available  to  be 
able  to  calculate  them  or  to  extract  the  Cl  energy  of  the  lowest  state.  As  a  result, 
one  is  often  forced  to  preselect  the  configurations  that  are  expected  to  contribute  the 
most  and  to  consider  only  those.  This  may  be  done  either  by  limiting  the  number 
of  electrons  excited  simultaneously  from  the  ground  state  HF  wave  function  or  by 
reducing  the  number  of  orbitals  from  which  or  into  which  excitations  are  permitted. 
The  limitations  imposed  are  then  included  when  referring  to  the  type  of  calculation 
performed. 

FULL-CI  considers  all  possible  excitations  from  the  HF  wave  function.  A 
FULL-CI  energy  represents  the  lowest  energy  one  can  attain  with  a  given  set  of 
orbitals.  More  often,  one  must  restrict  the  excitations  to  just  singles,  doubles,  triples, 
etc.,  or  combinations  of  these.  The  shorthand  technique  for  specifying  which  were 
considered  is  to  use  the  capital  letter  of  each  allowed.  Therefore,  a  CI-SDQ 
considers  all  possible  single,  double,  and  quadruple  excitations  from  the  HF 
configuration. 

Sometimes,  symmetry  can  be  used  to  reduce  the  number  of  orbitals  that  are 
treated  in  a  Cl  calculation.  For  example,  in  an  unsaturated  molecule,  one  might  only 
consider  excitations  between  orbitals  of  n  symmetry.  These  orbital  restrictions  are 
once  again  included  when  referring  to  the  calculations.  For  instance,  a-S,  jt-SD  Cl 
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specifies  only  single  excitations  from  molecular  orbitals  of  a  symmetry  and  single  and 
double  excitations  from  molecular  orbitals  of  n  symmetry. 

Section  1-3-2  M011er-Plesset  Perturbation  Theory 

A  non-variational  approach  to  electron  correlation  is  through  the  perturbative 
technique  developed  by  Miller  and  Plesset.12  They  express  the  effect  of  including 
electron  correlation  as  a  perturbation  to  the  Hamiltonian  derived  under 
Hartree-Fock  theory,  ft^.  Hence,  if  the  real  hamiltonian  is  Aq,  and  the  perturbative 
hamiltonian  is  ft,  with  a  weight  of  X,  then  they  are  related  by  the  following  equation 

Hq  -  Hhf  +  XH  (30) 

Under  Rayleigh-Schrodinger  perturbation  theory,  the  wave  function  then  becomes 

T0  -  V(hF  +  XY(1)  +  X2,P(2)  + ...  (31) 

and  the  resultant  energy  is13 

E0  -  E^l  +  XE(,)  +  X2Ea)  + ...  (32) 

where  is  called  the  ?h  order  correction  to  the  HF  wave  function,  is  the  ?h 
order  correction  to  the  HF  energy,  and  Y0  and  E^,  are,  respectively,  the  true  wave 
function  and  the  exact  energy. 
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An  important  practical  consideration  is  how  many  terms  are  included  in 
approximating  1^.  For  example,  an  MP2  calculation  includes  up  to  the  second  order 
term  in  the  M^ller-Plesset  perturbation  energy.  For  higher  levels,  such  as  fourth 
order,  not  all  of  the  lower  order  terms  are  always  considered.  For  instance,  the 
difference  between  MP4SDQ  and  MP4SDTQ  is  that,  in  the  former,  triple  excitations 
are  not  included,  so  that  expansion  is  actually  a  truncated  or  partial  fourth-order 
level  of  perturbation  theory. 

Perturbative  techniques  for  electron  correlation  provide  results  that  are 
size-consistent14,  whereas  Cl  techniques  less  than  full  do  not.  Perhaps  of  even 
greater  practical  importance,  perturbation  theory  requires  substantially  less  CPU 
time16  and  disk  space  than  do  Cl  techniques  that  include  the  same  levels  of 
excitation.  However,  the  drawback  to  perturbation  theory  is  that  molecules  which 
are  not  well  represented  by  a  single  Slater  determinant  are  poorly  described,  since 
the  perturbation  theory  expansion  is  slow  to  converge. 

Section  1-3-3  Multi-Configuration  SCF  Theory. 

Especially  when  attempting  to  calculate  diradical  species,  expressing  a  singlet 
wave  function  with  a  single  determinant  is  grossly  inadequate.  In  such  cases,  a 
multi-configuration  wave  function  must  be  used  as  a  zeroth  order  approximation. 
The  exact  number  of  configurations  required  to  describe  accurately  the  wave  function 
for  a  molecule  may  vary  from  state  to  state  and  also  for  a  given  state,  with  changes 
in  geometry. 
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In  cases  where  a  single-determinantal  SCF  wave  function  is  inadequate,  it  is 
usually  best  to  obtain  a  multi-configuration  (MC)  wave  function  by  performing  an 
MCSCF  calculation.  An  MCSCF  wave  function  has  the  same  form  as  a  Cl  wave 
function,  but  in  a  MCSCF  calculation,  both  the  orbitals  and  the  mixing  coefficients 
are  simultaneously  optimized.  Consequently,  for  the  same  set  of  configurations,  a 
MCSCF  calculation  provides  a  variationally  better  wave  function.  However,  Cl 
calculations  are  usually  capable  of  handling  many  more  configurations  than  MCSCF 
calculations. 

MCSCF  calculations  are  capable  of  providing  a  good  zeroth-order  wave 
function;  but,  for  all  but  atoms  and  very  small  molecules,  not  a  completely  correlated 
one.  Of  course,  additional  electron  correlation  can  be  added  to  improve  an  MCSCF 
wave  function,  just  as  correlation  can  be  added  to  improve  a  RHF  or  UHF  wave 
function.  In  the  event  that  additional  correlation  is  added,  excitations  from  all  the 
configurations  used  in  the  MCSCF  wave  function  are  usually  included. 

Section  1-4  Basis  Sets 

Of  equal  or  greater  importance  to  the  adequate  inclusion  of  electron 
correlation  is  the  adequate  description  of  the  orbitals  used  in  the  Slater  determinants 
that  comprise  SCF,  MCSCF  and  Cl  wave  functions.  If  MOs  are  written  as  LCAOs, 
the  quality  of  the  MOs  depends  on  the  quality  of  the  functions  used  to  represent  the 
atomic  orbitals.  These  functions  are  called  the  basis  set. 
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A  type  of  orbital  used  frequently  in  atomic  calculations  is  a  Slater  type 
orbitals.  These  orbitals  have  the  functional  form: 

<t>  =  Yta(e,o)  *  f(R,«)  *  e'CR  (33) 


where  is  the  angular  portion  of  the  wave  function,  f(R,£)  is  the  radial 

harmonic  portion,  and  R  is  the  radial  distance  from  atomic  center. 

The  drawback  of  Slater  type  orbitals  for  molecular  calculations  is  that 
multiplying  two  or  more  Slater  orbitals  on  different  centers,  gives  integrals  that  are 
difficult  to  evaluate.  Almost  all  calculations  today  are  done  using  basis  sets 
composed  of  several  gaussian  type  functions  in  place  of  the  exponential  part  of  Slater 
type  orbitals.  Mathematically, 

«-<*  -  (34) 

i-l 

where  aj  and  b;  are  constants  that  can  be  varied  to  give  the  closest  fit  to  a  Slater 
orbital  or  that  can  be  optimized  for  calculation  on  molecules. 

The  advantage  of  replacing  one  Slater  orbital  with  several  gaussians  is  that 
integrals  over  gaussians  can  be  easily  and  rapidly  evaluated,  since  the  product  of  two 
gaussians  on  different  centers  is  a  third  gaussian.16  However,  gaussians  are  not 
ideal  for  representing  Slater  orbitals,  since  Slater  orbitals  have  a  cusp  at  the  origin; 
whereas,  gaussians  do  not.  To  duplicate  an  entire  Slater  orbital  accurately  requires 
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many  gaussians.  Hence,  as  opposed  to  matching  the  entire  Slater  orbital,  many  basis 
sets  try  to  duplicate  the  portion  of  the  orbital  deemed  most  important.  For  inner 
shell  orbitals,  not  involved  much  in  bonding,  duplication  of  the  cusp  is  most 
important.  For  valence  shell  orbitals,  the  outer  portions  of  the  atomic  orbitals  are 
the  most  important,  since  they  are  involved  bonding. 

Examples  of  common  basis  sets  are  STO-3G,17  3-21G,18  6-31G*,19  DZ,20 
and  DZP.21  STO-3G  is  a  minimal  basis  set.  It  replaces  each  Slater  Type  Orbital 
(STO)  with  3  gaussians  (G).  A  minimal  basis  set  has  no  flexibility,  since  the  same 
gaussians  are  used  to  describe  all  orbitals  within  a  given  shell,  regardless  of  the 
environment  surrounding  the  atom  or  type  of  bond  that  an  orbital  is  used  to  form. 
Thus,  for  example,  the  2p  orbital  in  STO-3G  gives  a  much  better  description  of  C-C 
a  bonds  than  of  C-C  ir  bonds.22 

A  more  flexible  basis  set  is  3-2 1G.  It  is  a  split-valence  basis  set,  where  the 
inner  shell  orbitals  are  still  represented  by  3  gaussians,  but  the  valence  shell  orbitals 
are  represented  by  two  different  and  independent  groups  of  gaussians.  One  group 
has  2  gaussians,  optimized  for  short  distance  interactions;  the  other  consists  of  one 
gaussian,  optimized  for  longer  distance  interactions.  The  weights  of  the  two  groups 
of  gaussians  on  each  center  in  each  MO  are  determined  variationally  by  RHF,  UHF, 
or  MCSCF  calculations. 

The  6-31G*  basis  set  is  also  a  split-valence  basis  set.  Here,  however,  there 
are  6  gaussians  to  represent  inner  shell  orbitals  and  a  split  pair  of  3  and  1  gaussian 
functions  for  the  valence  orbitals.  The  asterisk  denotes  the  addition  of  polarization 
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functions  on  all  non-hydrogen  atoms.  This  mean  the  addition  of  d  orbitals  to  first 
and  second  row  atoms  and  the  addition  of  f  orbitals  to  atoms  below  row  two. 
Polarization  functions  allow  n  electron  density  to  be  shifted  away  from  the  center  of 
each  atom.  This  is  especially  important  in  small  rings  and  in  polar  bonds,  such  as  in 
carbonyl  groups,23  where  the  electron  density  is  concentrated  more  at  on  end  of  the 
bond  than  the  other. 

DZ  and  DZP  basis  sets  are  similar,  respectively,  to  3-21G  and  6-31G*,  but  in 
a  double  zeta  (DZ)  basis  set,  the  orbitals  in  each  shell,  not  just  the  valence  shell,  are 
all  represented  by  two  sets  of  gaussians.  The  most  commonly  used  DZ  and  polarized 
DZ  (DZP)  basis  sets  are  those  of  Dunning  and  Hay 20,21  They  found  that  different 
atoms  require  different  numbers  of  gaussians;  hence,  there  is  no  generic  specification 
of  the  number  of  gaussians,  as  there  is  in  the  basis  sets  derived  by  Pople  and 
coworkers  (e.g.  STO-3G,  3-21G,  and  6-31G*). 

Section  1-5  Isodesmic  Reactions 

Just  as  chemists  balance  equations  stoichiometrically,  computational  chemists 
often  balance  equations  isodesmically.24  By  ensuring  that  there  are  similar  types 
of  bonds  and  the  same  spin  multiplicity  of  states  on  both  sides  of  an  equation,  the 
magnitudes  of  the  errors  in  energies,  due  to  basis  set  deficiencies  and  inadequate 
treatment  of  electron  correlation,  can  be  reduced. 

For  example,  it  is  much  more  difficult  to  calculate  the  carbon-hydrogen  bond 
strergth  in  ethane  or  in  methane  exactly  than  to  calculate  the  carbon-hydrogen  bond 
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strength  in  ethane,  relative  to  that  in  methane.  Thus,  although  the  calculated  energy 
changes  for  (35)  and  (36)  may  not  accurate,  if  the  errors  in  each  are  identical,  then 
the  energy  for  (37)  will  be  exact,  since  this  isodesmic  reaction  is  simply  the  difference 
between  the  first  two. 


HCF^-CHj 

-»  H*  +  ‘CH^CH, 

(35) 

HCHj 

-  H*  +  *CH, 

(36) 

HCHj-CHj  +  ‘CH, 

-*  HCHg  +  ‘CHj-CHj 

(37) 

At  the  HF/6-31G*  level  of  theory,  the  calculated  bond  strength  of  the  C-H 
bond  in  ethane  is  73.4  kcal/mol  and  in  methane  is  76.1  kcal/mol.  These  values  are 
considerably  lower  than  the  experimental  values  of  100.4  kcal/mol,35  and  104.9 
kcal/mol26  respectively.  However,  since  equation  (37)  is  the  difference,  it  predicts 
that  the  C-H  bond  in  ethane  is  2.7  kcal/mol  weaker  than  in  methane,  which  is 
reasonably  close  to  the  experimental  difference  of  4.5  kcal/mol. 

The  way  to  improve  the  accuracy  of  the  energy  of  reaction  predicted  via  ab 
initio  techniques  is  to  increase  the  size  of  the  basis  set  or  increase  the  amount  of 
electron  correlation.  Inclusion  of  electron  correlation  through  second  order 
M011er-Plesset  perturbation  theory  usually  provides  a  closer  match  of  predicted  to 
experimental  results  than  do  calculations  at  the  SCF  level,  For  example,  using 
MP2/6-31G*  derived  energies,  the  C-H  bond  strength  in  ethane  is  90.7  kcal/mol  and 
in  methane  is  93.4  kcal/mol.  Inclusion  of  electron  correlation  at  the  MP2  level 
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obviously  provides  better  absolute  energies  of  reaction,  but  leaves  the  difference 
between  the  two  reactions  the  same,  2.7  kcal/mol.  MP4SDTQ/6-31G*  derived 
energies  give  only  a  small  further  increase  in  the  bond  strengths,  to  91.3  kcal/mol  for 
ethane  and  94.5  kcal/mol  for  methane.  The  energy  difference  also  increases  slightly 
to  3.1  kcal/mol.27 

The  absolute  bond  strengths  can  be  further  improved  through  addition  of 
more  extensive  electron  correlation  and  use  of  larger  basis  sets.28  However, 
isodesmic  reactions  allow  accurate  predictions,  usually  on  the  same  order  as  the 
experimental  uncertainty,  with  the  use  of  moderate  sized  basis  sets  and  low  levels  of 
electron  correlation.  Indeed,  in  the  example  discussed  above,  although  addition  of 
electron  correlation  dramatically  improves  the  calculations  for  the  absolute  C-H  bond 
strengths,  since  the  isodesmic  reactions  allows  for  cancellation  of  errors  in  each  of 
the  absolute  bond  strengths,  addition  of  electron  correlation  provides  no  or  little 
improvement  in  the  calculation  of  the  relative  C-H  bond  strengths. 

Isodesmic  reactions  are  useful  when  absolute  energies  are  of  less  interest  than 
differences  between  compounds.  In  those  cases  where  absolute  energies  are 
important,  they  can  be  calculated  via  isodesmic  reactions,  relative  to  an 
experimentally  known  value.  Because  calculations  of  absolute  values  for  quantities 
like  bond  dissociation  energies  are  computationally  demanding,  use  of  isodesmic 
reactions  provides  an  attractive  alternative. 
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Chapter  2  Substituent  Effects  on  Carbon-Hydrogen 
versus  Silicon-Hydrogen  Bond  Strengths 


Section  2-1  Experimental  Carbon-Hydrogen  versus  Silicon-Hydrogen  Bond  Strengths 
Carbon-centered  radicals  have  been  extensively  studied.1,2,3  Much  less  work 
has  been  done  on  silicon-centered  radicals.3,4,5  Since  carbon  and  silicon  are  both 
group  IVb  elements  in  the  periodic  table,  bonds  to  carbon  and  to  silicon  might  be 


presumed  to  be  similar.  However,  this  apparently  is  not  the  case. 


Table  2-1  Experimental  data*  for  bond  dissociation  energies  (kcal/mol)  for  the 
reaction,  R-H2X-H  -  R-H2X*  +  *H 

Compound 

X 

11 

O 

X  =  Si 

H3X-H 

104.8 

90.3 

(CH3)H2X-H 

100.4 

89.6 

CH3)2HX-H 

95.9 

89.6 

(CH3)3X-H 

93.3 

90.3 

(CH3)3C-H2X-H 

99.7 

(CH3)3Si-H2X-H 

99.2 

H3Si-H2X-H 

86.3 

[(CH3)3Si]3-X-H 

79.0 

“Data  from  reference  4,  except  for  the  BDE  of  CH3CH2-H,  which  comes  from 
reference  2,  and  [(CH3)3Si]3Si-H,  which  comes  from  reference  5. 
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Walsh4  found  experimentally  that  alkyl  substituents  have  different  effects  on 
silicon-hydrogen  versus  carbon-hydrogen  bond  strengths  (see  Table  2-1).  Alkyl 
substituents  appear  to  weaken  the  carbon-hydrogen  bond  by  roughly  4-5  kcal /mol  per 
substituent  group,  but  they  have  no  obvious  effect  on  silicon-hydrogen  bond  strengths. 
However,  silyl  substituents  weaken  both  silicon-hydrogen  and  carbon-hydrogen  bonds 
by  roughly  4  keal/mol  per  silyl  substituent. 

To  understand  these  differences  in  substituent  effects  on  silicon-hydrogen  and 
carbon-hydrogen  bonds  strengths,  we  decided  to  investigate  the  effects  of  a  wide 
range  of  substituents  by  computational  means.  We  examined  all  first  row 
substituents  plus  silicon,  phosphorous,  sulfur,  chlorine,  from  the  second  row  and  a 
limited  number  of  substituent  groups  containing  ir  bonds.  We  hoped  that  this  range 
of  substituents  would  be  sufficient  to  elucidate  completely  all  the  different 
mechanisms  by  which  substituents  affect  C-H  and  Si-H  bond  dissociation  energies 
(BDEs). 

Section  2-2  Theoretical  Methodology 

To  calculate  substituent  effects  on  carbon-hydrogen  and  silicon-hydrogen  bond 
strengths,  we  used  equation  2.1.  Equation  2.1  is  an  isodesmic  reaction  that  measures 


X-AHj  +  AH,  -*  X-AH,  +  AH, 


(2.1) 
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the  effect  of  substituent  group  X  on  the  A-H  bond  strength,  relative  to  that  in 
methane  (A  =  C)  or  silane  (A=Si).  As  discussed  in  Chapter  1,  use  of  an  isodesmic 
reaction  allows  cancellation  of  most  of  the  errors  made  in  calculation  of  the 
individual  BDEs, 

The  radical  geometries  were  optimized  with  a  UHF  wave  function  and  the 
6-31G*  basis  set.  Initially,  the  geometries  were  fully  optimized  with  respect  to  all 
degrees  of  freedom.  Later,  some  radicals  were  also  optimized  to  a  transition  state 
involving  a  conformational  change.  Many  of  these  transition  states  were  located  by 
imposing  symmetry  constraints.  Others  were  optimized  with  the  constraint  that  the 
second  derivative  of  the  Hessian  matrix  have  exactly  one  negative  eigenvalue. 
Electron  correlation  was  then  added  using  MP4SDTQ  at  the  UHF  optimized 
geometry.  The  MP4SDTQ/6-31G*  energies  for  the  reaction  of  equation  (2.1)  are 
given  in  Tables  2-2  through  2-6. 

The  UHF  wave  functions  for  the  allylic  type  radicals  proved  to  be  have  a  large 
degree  of  spin  contamination.  In  order  to  obtain  a  pure  doublet  wave  function, 
instead  of  using  MP4SDTQ  calculations  to  include  electron  correlation,  we  employed 
the  direct  CI-SEP  routine  in  GAUSSIAN  867  This  CI-SD  routine  calculates  the  Cl 
wave  function  energy  considering  all  single  and  double  excitations  from  the  lowest 
energy  configuration.  The  CI-SD  energies  for  the  allylic  systems  are  given  in  Tables 
2-7  through  2-9. 
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Section  2-3  First  Row  Substituents 

First  row  substituents  can  be  classified  into  three  distinct  types:  (1)  those  with 
empty  p  orbitals,  (2)  those  with  lone  pairs  of  electrons,  and  (3)  those  without  either, 
e.g.  methyl.  Each  of  these  groups  of  substituents  behave  differently  and  so  will  be 
discussed  individually. 

Section  2-3-1  Substituents  with  Empty  p  Orbitals 

The  early  first  row  substituents,  Li,  BeH  and  BHg,  all  have  at  least  one  empty 
p  orbital  available  to  delocalize  the  odd  electron  at  a  carbon  or  silicon  radical  center. 
The  resulting  partial  n  bonds  should  stabilize  the  substituted  silyl  or  methyl  radicals. 
Indeed,  as  shown  in  Table  2-2,  equation  2.1  is  calculated  to  be  exothermic  in  all  cases 
for  both  carbon  and  silicon.  This  indicates  that  all  of  these  substituents  lower  the 
bond  dissociation  energy  (BDE)  of  the  A-H  bond. 

Testing  for  the  existence  of  n  delocalization  in  X-Af^  as  teh  probable  source 
of  stabilization  is  straightforward  using  the  BHj  substituent.  Rotation  of  the  empty 
p  orbital  out  of  conjugation  with  the  singly  occupied  orbital  on  the  radical  center 
precludes  any  stabilization  through  n  delocalization.  The  barrier  to  rotation  about 
the  carbon-boron  bond  in  HgC-BHj  is  9.7  kcal/mol.  Since  the  boron  substituent 
reduces  the  C-H  BDE  by  12.4  kcal/mol,  most  of  the  reduction  in  the  BDE  by  BE^ 
comes  from  delocalization  of  the  odd  electron  on  carbon  into  the  empty  p  orbital  on 


boron. 


Table  2-2  Effect  of  substituents  with  empty  p  orbitals  on  the 
MP4SDTQ/6-31G*  derived  BDEs  using  the  reaction*  [X-AH^*  +  AH*  -» 
X-AH,  +  .AH, 


X-A  bond  length  in  the 
AE  (kcal/mol)  radical  (A) 


Substituent 
Group  (X) 

u 

II 

< 

c/3 

n 

< 

A  =  C 

A  =  Si 

Li 

9.4 

12.0 

1.943 

2.426 

BeH 

9.0 

8.0 

1.665 

2.166 

BHj15 

12.4 

12.2 

1.536 

1.954 

BH/ 

2.7 

0.8 

1.555 

2.046 

BlV 

2.7 

-4.5 

1.555 

2.045 

aA  positive  AE  implies  that  the  substituent  weakens  the  A-H  bond. 
bBoth  carbon  and  silicon  are  planar. 

'Radical  center  is  rotated  out  of  conjugation  with  empty  p  orbital  and  allowed 
to  pyramidalize.  Carbon  remains  planar,  but  silicon  has  an  out  of  plane  angle 
(0)  of  51°. 

d Radical  center  is  rotated  out  of  conjugation  with  empty  p  orbital  and 
constrained  to  be  planar. 

• 

The  BjB-SiHj  radical  gives  similar  results.  Rotation  about  the  B-Si  bond 
requires  11.4  kcal/mol.  Upon  rotation  however,  there  is  significant  change  in  the 
geometry  at  the  silicon  center.  In  the  conjugated  geometry,  the  silicon  center  is 
planar;  but  upon  rotation,  the  silicon  becomes  significantly  pyramidal.  Therefore,  the 
energy  change  is  not  due  strictly  to  the  breaking  of  the  n  bond  but  also  contains  the 
energy  lowering  that  accompanies  pyramidalization  of  the  silicon  center.  Hence,  to 
get  a  true  measure  of  the  silicon-boron  ir  delocalization  energy,  the  silicon  center 
must  be  held  planar  when  rotated.  Optimization  of  this  Qv  structure  reveals  that 
breaking  the  ir  bond  requires  16.8  kcal/mol,  and  pyramidalization  regains  5.4 
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kcal/mol  of  that  energy.  The  larger  magnitude  of  the  ir  delocalization  energy  in 

•  • 

HjB-SiH;  than  in  HjB-CHg  effect  can  be  understood  using  either  resonance  theory 
or  qualitative  MO  perturbation  theory. 

In  resonance  theory,  delocalization  of  the  odd  electron  on  A  by  an  empty 
orbital  on  X  is  represented  using  equation  2.2.  The  amount  of  stabilization  should 


X-AHj  —  X-A+H,  (2.2) 

be  proportional  to  the  contribution  the  second  resonance  structure  makes.  The 
second  resonance  structure  should  contribute  mre  as  the  electronegativity  of  atom 
A  decreases.  This  explains  why,  although  silicon  generally  forms  weaker  n  bonds 
than  carbon,8  the  n  bond  strength  is  greater  in  H^B-AL^  for  A  =  Si  than  for  A  =  C. 

The  second  resonance  structure  should  also  contribute  more  as  the 
electronegativity  of  X  increases.  Hence,  the  anticipated  order  of  ir  bond  strengths 
for  substituent  groups,  X,  is:  BHj  >  BeH  >  Li.  However,  this  is  not  the  calculated 
order  found  in  Table  2.2.  The  actual  order  has  Li  lowering  the  A-H  BDE  more  than 
BeH  for  both  A=C  and  A  =  Si.  Apparently,  there  is  a  second  effect  that,  unlike  the 
formation  of  the  1-electron  ir  bond,  increases  the  A-H  BDE  as  the  electronegativity 
of  the  substituent  group,  X,  increases. 

Use  of  MO  perturbation  theory  gives  an  identical  prediction  of  the  effect  of 
substituent  electronegativity  on  the  strength  of  the  1-electron  n  bond  in  X-AH,. 
Since  the  empty  orbital  on  X  and  the  orbital  containing  the  odd  electron  on  A 
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overlap,  there  will  be  mixing  of  these  two  orbitals.  Using  second-order  perturbation 
theory,  the  amount  of  energy  lowering  that  results  is  roughly  proportional  to  the 
square  of  the  overlap  of  these  two  orbitals,  divided  by  the  energy  gap  between  them. 
The  less  electronegative  A  is,  the  closer  the  half  filled  orbital  is  to  the  empty  orbital 
in  energy,  the  smaller  the  energy  gap  between  them,  and  the  larger  the  energy 
lowering  due  to  their  mixing.  Once  again,  this  explains  why  the  n  bond  to  BHj  is 
stronger  for  A=Si  than  A  =  C. 

The  more  electronegative  the  substituent  group  X  is,  the  closer  in  energy  its 
empty  p  orbital  is  to  the  half  filled  orbital,  and  once  again,  the  larger  the  energy 
lowering  caused  by  the  formation  of  the  resultant  ir  bond.  MO  perturbation  theory 
thus  also  predicts  that  the  anticipated  order  for  the  substituent  groups  in  decreasing 
the  A-H  BDE  is:  X  =  BF^  >  X  =  BeH  >  X=Li.  As  indicated  in  the  paragraph  above, 
the  data  in  table  2.2  shows  that  there  is  an  additional  effect,  not  yet  identified,  that 
apparently  acts  to  increase  the  A-H  BDE  as  the  electronegativity  of  X  increases. 

Section  2-3-2  Substituents  with  Lone  Pairs 

The  later  first  row  substituents,  X  =  NHj,  OH,  and  F,  all  contain  at  least  one 
lone  pair  of  electrons  that  can  interact  with  the  radical  center  on  A.  Stabilization  of 
the  odd  electron  on  A  can  occur  through  formation  of  a  3-electron  jt  bond. 
Although  one  of  the  three  electrons  will  be  in  an  anti-bonding  n  orbital,  the  two 
electrons  that  occupy  the  bonding  *  orbital  should  result  in  net  stabilization  for  the 
3-electron  system.  Since  Nf^  has  a  single  lone  pair,  the  strength  of  the  3-electron 
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Table  2-3  Effect  of  substituents  with  lone  pairs  on  the  MP4SDTQ/6-31G* 
derived  BDEs  using  the  reactioif  [X-AP^]*  +  AH*  -»  X-AT^  +  *AHg 

AE  (kcal/mol) 

X-A  bond  length  in  the 
radical  (A) 

Substituent 

A  P 

A  =  Si 

a  r 

A  =  Si 

Group  (X) 

A-L 

Pi.  — 

NHj 

12.2 

0.9 

1.40? 

1.7291 

NH* 

3.7 

-2.5 

1.420s 

1.721f 

OH 

8.6 

-0.8 

1.359 

1.652 

OH15 

4.2 

-1.6 

1.366 

1.652 

F 

4.4 

-3.5 

1.331 

1.599 

“A  positive  AE  implies  that  the  substituent  weakens  the  A-H  bond. 
bRadical  center  is  rotated  out  of  conjugation  with  lone  pair  orbital. 
'Pyramidalization  angle  at  C  (< )  is  34*,  =46°. 

d<fci  =  56*,0*  =  28*. 

*4c  =  4.4‘,<Ps  =  52°. 
f<fci  =  48°,  (pfi  =  10*. 


7r  bond  in  HjN-AHj  can  be  measured  by  rotation  of  the  substituent  group  out  of 
conjugation  with  the  radical  orbital. 

Our  calculations  find  that  in  RjN-CHj,  9.5  kcal/mol  is  required  to  rotate  the 
carbon  radical  out  of  conjugation  with  the  nitrogen  lone  pair.  Even  at  the  rotated 
geometry,  hyperconjugation  of  the  radical  with  the  N-H  bonds  of  the  amino  group 
may  provide  some  stabilization  for  the  radical.  Thus,  the  9.5  kcal/mol  barrier  to 
NHj  rotation  indicates  that  most,  if  not  all,  of  the  12.2  kcal/mol  reduction  in  the  C-H 
BDE  by  NHj  is  due  to  formation  of  a  n  bond  in  the  radical. 

Additional  evidence  of  the  formation  of  a  partial  w  bond  in  HjN-CHj  comes 
from  careful  examination  of  the  geometry  changes  that  occur  upon  rotation  of  the 
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I^N  group.  The  complete  geometries  are  given  in  Appendix  I,  but  some  parameters 
are  listed  in  Table  2-3.  In  the  equilibrium  structure,  the  N-C  bond  length  is  1.402A, 
which  is  shorter  than  the  N-C  bond  length  of  1.420A  in  the  rotated  structure. 

At  first  glance,  since  neither  atomic  center  is  planar  in  either  of  these 
structures,  it  is  not  obvious  that  a  rr  bond  that  is  formed;  but  the  changes  in 
pyramidalization  at  the  two  centers  are  significant.  In  the  equilibrium  geometry,  the 
pyramidalization  angles  are  34°  for  carbon  and  46°  for  nitrogen;  for  the  rotated 
structure,  they  are  4.4  °  for  carbon  and  52  °  for  nitrogen.  The  carbon  center  is  more 
pyramidal  and  the  nitrogen  is  less  pyramidal  in  the  equilibrium  structure  because  of 
a  transfer  in  electron  density  from  the  nitrogen  to  the  carbon,  as  depicted  in  the 
second  resonance  structure  of  equation  2.3.  An  increase  in  electron  density  at  the 

H^N -CR,  «_  (2.3) 

carbon  center  would  tend  to  make  the  carbon  center  behave  more  like  a  carbanion, 
and  carbanions  generally  prefer  pyramidal  geometries.  The  reduction  in  electron 
density  at  the  nitrogen  would  tend  to  flatten  it.  Hence,  the  changes  in  geometry  are 
consistent  with  formation  of  a  3-electron  n  bond. 
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Both  resonance  theory  and  MO  perturbation  theory  predict  that  formation  of 
a  three-electron  tt  bond  requires  electron  transfer  from  the  substituent  group,  X,  to 
A.  In  resonance  theory,  the  transfer  is  represented  in  the  second  resonance  structure 


X-AHj  — *  X+-A  Hj  (2.4) 

in  equation  2.4.  In  MO  theory,  when  the  orbitals  on  X  and  A  mix,  some  of  the 
electron  density  in  the  doubly  occupied  orbital  on  X  is  transferred  to  the  singly 
occupied  orbital  on  A.  Using  either  of  these  theories,  the  greatest  stabilization  of 
the  radical  and  hence,  the  greatest  lowering  of  the  BDE  is  expected  when  the 
difference  in  electronegativity  between  X  and  A  is  smallest.  This  predicts  that 
carbon,  which  is  more  electronegative  than  silicon,  should  have  a  stronger  n  bond  in 
X-AHj  and  hence  a  lower  A-H  BDE  than  silicon.  Also,  the  anticipated  trend  for  the 
n  bond  strengths  of  the  substituent  groups,  X,  should  be  NHj  >  OH  >  F.  This  is 
exactly  what  is  found  in  the  calculated  A-H  BDEs  that  are  listed  Table  2-3. 

However,  the  Si-H  BDE  in  X-AHg  actually  increases  relative  to  that  in  silane 
for  X  =  OH,  and  F,  and  for  X  =  NHj  when  the  amino  lone  pair  is  rotated  out  of 
conjugation.  Three  electron  n  bonding  in  the  radical  obviously  cannot  explain  these 
increases  in  BDE.  Hence,  once  again,  there  is  evidence  of  another  effect  acting;  and 
once  again,  this  effect  appears  to  increase  the  Si-H  BDEs  as  the  electronegativity  of 
X  increases.  Thus,  while  NF^,  OH  and  F  can  stabilize  silicon  centered  radicals 
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through  ir  electron  donation,  these  substituents  increase  the  Si-H  BDE  presumably 
by  withdrawing  electron  density  from  silicon  through  the  a  bond.  There  is 
experimental  evidence  for  the  existence  of  this  second  effect  in  both  Si-H  and  C-H 
BDEs. 

Experimentally,  it  is  found  that  in  FSA-H,  the  A-H  bond  is  stronger  than  in 
AH,  for  both  A=C  and  A  =  Si.  The  Si-H  bond  strength  in  trifluorosilane  has  been 
measured  at  100.1  kcal/mof  which  is  9.8  kcal/mol  greater  than  that  in  silane.  In 
addition,  despite  the  probable  stabilization  of  CFS  The  C-H  bond  strength  of  106.0 
kcal/mol10  in  trifluoromethane  is  greater  than  that  in  methane  by  1.2  kcal/mol. 

Additional  evidence  exists  that  indicates  C-H  bond  strengths  can  be  increased 
by  a  single  substituent  group  that  is  sufficiently  electron  withdrawing.  Thus, 
experimentally,  1,1,1-trifluoroethane  has  a  carbon-hydrogen  bond  strength  of  106.2 
kcal/mol,1  presumably  because  of  the  presence  of  the  very  electron  withdrawing  CFS 
group  attached  to  the  incipient  carbon  radical  center.  We  find  that  the  C-H  BDE 
in  HgN*  -CHg  exceeds  that  in  CH<  by  4.0  kcal/mol,  again  because  the  H  N+  group 
is  very  electron  withdrawing. 

Although  both  C-H  and  Si-H  BDE  are  increased  by  a  electron  withdrawal, 
this  effect  is  masked  by  n  delocalization  in  :X-CH-  Since  n  delocalization  is  less 
important  in  :X-SiH»  the  a  withdrawing  nature  of  substituents  more  electronegative 

than  hydrogen  actually  results  in  an  increase  in  calculated  BDEs  for  :X-SiHg,  except 

• 

for  X  =  NHj  when  the  nitrogen  lone  pair  is  allowed  to  conjugate  in  the  F^N-SiHj 


radical. 
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The  a  electronegativity  effect  also  explains  why  X  =  Li  reduces  both  C-H  and 
Si-H  BDEs  more  than  X  =  BeH.  Lithium  is  more  electropositive  than  beryllium;  and, 
hence,  although  Be  is  a  better  n  acceptor,  Li  is  a  better  a  electron  donor.  It  appears 
that  the  a  electron  donation  effect  dominates  for  these  two  substituents. 

For  the  early  first  row  substituents,  BF^  is  actually  an  aberration.  Its 
electronegativity  allows  it  to  be  very  effective  at  delocalizing  the  odd  electron  on  A 
by  forming  a  strong  n  bond  to  both  A  =  C  and  Si.  This  masks  the  fact  that  it  is  a 
worse  a  electron  donor  than  either  BeH  or  Li. 

Thus,  our  calculations  on  n  acceptors  and  *  donor  substituents  reveal  that 
there  are  actually  two  types  of  substituent  effects  on  A-H  BDE.  One  is  electron 
delocalization  by  ir  bonding  in  the  radical,  which  always  lowers  the  A-H  BDE.  The 
other  is  a  a  effect,  which  appears  to  lower  the  A-H  BDE  when  the  substituent,  X, 
is  less  electronegative  than  hydrogen  and  raises  the  BDE  when  X  is  more 
electronegative  than  hydrogen. 

Section  2-3-3  Methyl  and  Silyl  as  Substituents 

Based  on  these  two  effects,  we  are  now  prepared  to  answer  the  original 
question:  Why  do  silyl  groups  weaken  both  silicon-hydrogen  and  carbon-hydrogen 
bonds  whereas  alkyl  groups  only  weaken  carbon-hydrogen  bonds? 

Granted,  n  delocalization  is  less  important  for  substituents  like  methyl  and 
silyl  that  do  not  have  a  lone  pair  or  an  empty  p  orbital;  but  hyperconjugation  does 
allow  delocalization  of  the  odd  electron,  as  depicted  in  equations  2.5  and  2.6. 
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Table  2-4  Effect  of  methyl  and  silyl  substituents  on  the  MP4SDTQ/6-31G* 
derived  BDEs  using  the  reaction^  [X-AHj]»  +  AH<  -»  X-AH^  +  »Ab^ 

Substituent 
Group  (X) 

AE  (kcal/mol) 

X-A  bond  length  in  the 
radical  (A) 

A  =  C 

> 

ii 

<z> 

U 

ii 

< 

A=Si 

CI^ 

3.3 

-0.8 

1.498 

1.894 

SiHj 

4.4 

2.9 

1.860 

2.345 

*A  positive  AE  implies  that  the  substituent  weakens  the  A-H  bond. 


Hyperconjugation  should  play  a  larger  role  in  stabilizing  carbon-centered  than 
silicon-centered  radicals,  since,  the  resultant  ir  bond  to  carbon  is  stronger  than  the 
k  bond  to  silicon  for  both  X  =  Cf^  and  X  =  SiH^.  Thus,  hyperconjugation  is  expected 
to  have  a  larger  effect  on  reducing  the  C-H  BDE  than  the  Si-H  BDE. 

It  is  a  little  harder  to  predict  whether  a  substituent  should  hyper  conjugate 
better  than  a  SiH^,  because  the  difference  in  the  strengths  of  the  *  bonds  formed  is 
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partially  compensated  for  by  the  difference  in  the  strengths  of  the  bond  that  is 
broken.11  However,  because  SiHg,  unlike  CH,,  is  less  electronegative  than 
hydrogen,  the  a  donating  effect  of  X  =  SiHg  should  lower  the  BDE  for  both  A  = 
C  and  Si. 

As  expected,  the  results  in  Table  2-4  show  that  both  substituents  reduce  the 
C-H  BDE.  Despite  the  fact  that  CHg  hyperconjugation  is  expected  to  be  slightly 
more  radical  stabilizing  than  SiHg  hyperconjugation,11  an  SiHg  substituent  is 
calculated  to  cause  a  slightly  larger  reduction  in  the  C-H  BDE.  Presumably  this  is 
due  to  the  fact  that  SiHg  is  a  better  a  donor  than  CHg. 

Because  of  the  weakness  of  n  bonds  to  silicon,  neither  substituent  is  expected 
to  reduce  the  Si-H  BDE  significantly  through  hyperconjugative  radical  stabilization. 
Consequently,  since  methyl  is  slightly  more  a  electron  withdrawing  than  hydrogen, 
it  is  perhaps  not  entirely  surprising  that  a  methyl  substituent  is  actually  calculated  to 
slightly  increase  the  Si-H  BDE.  In  contrast,  the  a  donating  ability  of  the  SiHg 
substituent,  relative  to  H,  results  in  the  lowering  of  the  Si-H  BDE  in  HgSi-SiH, 
relative  to  that  in  SiH,. 

Section  2-3-4  Origin  of  the  a  Electronegativity  Effect 

In  addition  to  the  well  understood  effect  of  substituent  electronegativity  on  jr 
delocalization  in  X-AHj  radicals,  our  calculations  have  revealed  that  a  electron 
donor  substituents  weaken  A-H  bonds  while  a  electron  withdrawing  substituents 
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strengthen  them.  This  second  effect  has  been  little  appreciated,  but  it  can  be 

understood  on  the  basis  of  a  rule  enunciated  by  Bent: 

Atomic  p  character  concentrates  in  orbitals  directed  toward  electronegative 
substituents.12 

Concentration  of  the  p  character  of  atom  A  in  the  X-A  bond  in  X-AHg  implies  that 
s  character  will  be  concentrated  in  the  A-H  bonds.  For  an  electronegative 
substituent,  X,  this  should  strengthen  the  A-H  bonds,  since  the  bond  strength 
increases  with  increasing  s  character  on  A.13  Conversely,  when  X  is  more 
electropositive  than  H  s  character  will  be  concentrated  in  the  X-A  bond. 
Consequently,  p  character  will  be  concentrated  in  the  A-H  bonds  and  weaken  them. 

The  a  electronegativity  effect  is  particularly  apparent  in  Si-H  BDE  when 
electronegative  substituents  are  present  since,  as  discussed  above,  ir  delocalization 
by  donation  of  substituent  lone  pairs  is  not  very  important  for  silicon-centered 
radicals.  Similarly,  since  hyperconjugation  provides  little  stabilization  for 
silicon-centered  radicals,  the  a  electronegativity  effect  is  particularly  apparent  in  the 
comparison  of  how  alkyl  and  silyl  substituents  affect  Si-H  BDEs. 

Section  2-3-5  First-row  Substituents  Containing  it  Bonds 

Organic  chemists  have  long  known  that  adjacent  n  bonds  stabilize 
carbon-centered  radicals  by  n  electron  delocalization14.  From  the  preceding 
discussion,  adjacent  n  bonds  should  stabilize  carbon-centered  more  than 
silicon-centered  radicals,  since  silicon  forms  weaker  n  bonds.  To  test  this  hypothesis, 
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we  examined  computationally  the  effect  of  several  different  substituent  groups  that 
are  commonly  used  to  stabilized  carbon  centered  radicals. 


Table  2-5  Effect  of  first  row  substituents  containing  n  bonds  on  the 
MP4SDTQ/6-31G*  derived  BDEs  using  the  reactiorf  [X-AI-^]*  +  AH<  — * 
X-AHg  +  .AH, 


AE  (kcal/mol) 


Substituent 

Group  (X) 

A  =  C 

A  =  Si 

-CN 

5.1 

-1.0 

-NC 

7.4 

-1.6 

-QH3 

13.2 

-3.1 

-CHO 

5.4 

0.4 

-NO, 

2.5 

-1.2 

aA  positive  AE  implies  that  the  substituent  weakens  the  A-H  bond. 


As  expected,  all  the  substituents  in  Table  2-5  are  calculated  to  lower  the  C-H 
BDE.  None  of  these  substituents  reduced  the  Si-H  BDE  from  that  in  SiH,,  since  ir 
delocalisation  is  relatively  ineffective  in  stabilizing  silyl  radicals;  and  all  of  these 
substituents  are  relatively  a  electron  withdrawing.  The  effect  of  substituents 
containing  n  bonds  to  second-row  elements  will  be  discussed  in  Section  2-5. 

Section  2-4  Second  Row  jt  Donor  Substituents 

Second  row  elements  differ  from  their  corresponding  first  row  elements  by 
forming  weaker  n  bonds  and  by  being  more  electropositive.  Since  Si-H  BDEs  are 
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more  affected  by  the  a  electronegativity  of  substituents  than  by  their  bonding 
abilities,  second  row  substituents  should  reduce  Si-H  BDEs  more  than  their  first  row 
counterparts. 


Table  2-6  Comparison  of  effects  of  second  and  first  row  substituents  on  the 
MP4SDTQ/6-31G*  derived  BDEs  using  the  reactioif  [X-AF^*]  +  AE^  -♦ 
X-AHg  +  *AHg 


AE  (kcal/mol)  AE  (kcal/mol) 


Subst. 
Group  (X) 

> 

it 

n 

> 

n 

C/5 

Subst. 
Group  (X) 

A=C 

> 

n 

C/5 

PH, 

5.7 

2.6 

NH, 

12.2 

0.9 

PH,b 

2.8 

0.7 

NH," 

3.7 

-2.5 

PH,' 

10.4 

NHj' 

12.8 

SH 

7.9 

1.6 

OH 

8.6 

-0.8 

SH* 

3.6 

-0.2 

OH15 

4.2 

-1.6 

Cl 

4.9 

-0.3 

F 

4.4 

-3.5 

aA  positive  AE  implies  that  the  substituent  weakens  the  A-H  bond. 
bRadical  orbital  rotated  out  of  conjugation  with  lone  pair. 
'Substituent  group  constrained  to  be  planar. 


The  expected  effect  of  these  substituents  on  the  C-H  BDE  is  more 
complicated.  The  lower  electronegativity  of  second  rows  substituents  should  make 
them  better  n  donors  than  their  first  row  counterparts,  but  the  weaker  n  bonds  that 
second  row  substituents  form  could  cause  them  to  provide  less  stabilization  for  the 
radicals  than  their  first  row  counterparts,  hence  increasing  the  C-H  BDE.  Of  course, 
the  lower  electronegativity  of  the  second  row  substituents  should  again  result  in  the 
a  effect  acting  to  lower  the  C-H  BDEs,  relative  to  their  first  row  counterparts. 
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The  second  row  substituents  on  silicon  behave  just  as  expected.  Each  gives 
a  lower  calculated  Si-H  BDE  than  its  first  row  counterpart.  X  =  PH2  lowers  the  Si-H 
BDE  more  than  X  =  NE^,  despite  the  fact  that  the  latter  substituent  provides  more 
■n  bonding  in  the  X-SiHj  radical.  Evidence  for  a  stronger  it  bond  in  X  =  Nf^  comes 
from  the  finding  that  it  costs  1.9  kcal/mol  to  rotate  the  group  out  of  conjugation, 
compared  to  3.6  kcal/mol  for  NHj. 

For  carbon,  SH  and  Cl  substituents  are  calculated  to  give  almost  the  same 
reduction  in  C-H  BDE  as  their  first  row  counterparts.  PHj  however  provides  much 
smaller  reduction  than  NF^. 

The  much  smaller  barrier  to  rotating  the  PF^  out  of  conjugation  in  X-CF^ 
again  suggests  that  ir  delocalization  in  the  radical  is  much  less  stabilizing  for  PHj 
than  NFL,.  This  is  not  surprising  since  the  phosphorous  lone  pair  does  not  overlap 
well  with  the  carbon  radical  orbital,  due  to  the  highly  pyramidal  geometry  at 

phosphorous,  which  results  in  a  large  amount  of  s  character  in  the  phosphorous  lone 

• 

pair.  If  the  phosphorous  is  constrained  to  be  planar  in  both  H^P- AHg  and  HjP-AF^, 
the  amount  of  it  stabilization  increases.  This  is  evidenced  by  the  fact  that  a  planar 
PFLj  substituent  is  calculated  to  reduce  the  C-H  BDE  by  10.4  kcal/mol.  The  C-H 
BDE  reduction  and  barrier  to  rotation  for  a  planar  NHj  substituent  are  12.8  and  10.3 
kcal/mol. 

In  contrast  to  PF^,  both  SH  and  Cl  substituents  have  lone  pairs  in  pure  p 
orbitals.  This  accounts  for  the  fact  that  they  reduce  the  C-H  BDE  more  than  a 
pyramidalized  PF^  group  but  not  as  much  as  a  planar  PHj.  The  barrier  of  4.3 
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kcal/mol  to  rotation  of  the  lone  pair  in  the  pure  p  orbital  on  S-H  out  of  conjugation 
is  also  intermediate  between  the  values  of  2.9  and  21.2  kcal/mol  for  rotating, 
respectively,  pyramidalized  and  planar  PH,  groups  in  PH^-CH,. 

Section  2-5  Vinylic  Substituents  Containing  Silicon 

Organic  chemists  have  long  known  that  vinyl  substituted  carbon  radicals  and 
ions  possess  additional  stability  compared  to  alkyl  substituted  carbon  radicals  and 
ions.1  This  additional  stability  comes  from  the  ability  of  the  n  bond  in  vinyl  to 
delocalize  the  radical  or  charge  through  allylic  resonance,  as  shown  in  equation  2.7. 
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(2.7) 


We  wanted  to  know  if  allylic  type  resonance  would  bestow  comparable 
stability  and  delocalization  when  silicon  replaces  one  or  more  of  the  carbon  atoms. 
Therefore,  we  performed  calculations  on  the  effect  of  vinyl  substituents  and  their 
H^C=SiH,  H,Si=CH,  and  H,Si  =  SiH  analogs  on  C-H  and  Si-H  BDEs. 
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Since,  as  discussed  above,  substituents  alter  BDE  through  both  a  and  n 
effects,  we  separated  these  effects  by  the  now  familiar  technique  of  rotating  the  vinyl 
group  out  of  conjugation  in  the  radical. 

Section  2-5-1  Computational  Methodology  Used 

A  problem  with  UHF  calculations  is  that  they  give  wave  functions  for  radicals 
that  are  not  pure  doublet  states.  In  X-AHj  radicals,  the  spin  contamination  becomes 
serious  when  X  is  unsaturated.  For  example,  for  allyl  radical  (X  =  I^C=CH  and 
A  =  C)  a  UHF  calculation  with  the  6-3 1G*  basis  set  gives  a  wave  function  with  J? 
=  0.972,  compared  to  S?  =  0.750  for  a  pure  doublet. 

When  electron  correlation  is  provided  by  M^ller-Plesset  perturbation  theory, 
spin  contamination  raises  the  energy  of  the  resulting  wave  function,  since  the 
contaminant  states,  quartets  and  even  higher  multiplicities,  lie  above  the  doublet 
wave  function  in  energy.  Since  spin  contamination  is  much  more  of  a  problem  in 
X-AHj  radicals  than  in  AHg,  especially  when  X  is  unsaturated,  the  values  for  the 
lowering  of  the  A-H  BDEs  by  the  unsaturated  substituents  in  Table  2-5  should  be 
regarded  as  lower  bounds  on  the  actual  BDEs. 

One  way  of  handling  spin  contamination  is  to  project  the  largest  spin 
contaminant  from  a  UHF  wave  function.15  However,  rather  than  using  spin 
projection,  we  elected  to  perform  CI-SD  calculations,  which,  being  strictly  variational, 
necessarily  give  pure  doublet  wave  functions.  The  effects  of  the  vinyl  substituents  on 
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the  A-H  BDEs  in  Tables  2-7  and  the  barriers  to  rotation  in  Table  2-8  were  all 
calculated  from  CI-SD  energies. 

Section  2-5-2  Effects  of  Vinyl  Substituents  Containing  Silicon  on  BDEs 

Since,  as  discussed  above,  n  bonding  to  C  is  generally  more  important  than 

it  bonding  to  Si,  one  would  expect  that  each  type  of  vinyl  substituent,  X,  would  have 

•  • 

a  greater  stabilizing  effect  in  X-CHj  than  in  X-SiH^  Therefore,  one  would  expect 
that  in  all  cases,  the  C-H  bind  in  X-CHj  would  be  weakened  significantly  more  than 
the  Si-H  bond  in  X-SiHj. 

Table  2-7  Effects  of  vinyl  derivatives  on  the  MP4SDTQ/6-31G*  derived  BDEs 
using  the  reaction*  [X-AHJ*  +  AH<  -♦  X-AH,  +  »AH, 


Vinyl 

Substituent 

Group  (X) 

AE  (kcal/mol) 

A=C 

A=Si 

I^C=CH 

17.2 

-0.7 

I^Si=CH 

36.3 

-0.2 

E^C=SiH 

13.1 

3.3 

I^Si  =  SiH 

18.8 

8.4 

*A  positive  AE  implies  that  the  substituent  weakens  the  A-H  bond. 

As  shown  in  Table  2-7,  this  is  the  case.  However,  the  amount  of  reduction  in 
BDE  reduction  depends  on  the  identity  of  the  vinyl  group,  X,  in  X-AHg,  as  well  as 
on  the  identity  of  A  Interestingly,  HjSi=CH  has  by  far  the  greatest  effect  of  all  the 
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substituents  on  reducing  the  C-H  BDE,  but  X  =  E^Si  =  SiH  has  the  largest  effect  on 
reducing  the  Si-H  BDE. 

Since  it  delocalization  in  the  allylic  radicals  formed  is  anticipated  to  be  the 
dominant  effect  on  C-H  BDEs,  it  should  be  possible  to  understand  the  effects  of  the 
various  vinyl  groups  by  examining  the  two  allylic  resonance  structures  in  equation  2.8. 
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For  Y=C,  Z=C  or  Si,  and  A=C,  the  two  structures  are  identical  and,  hence, 
equally  important.  However,  for  Y=Si  and  Z  =  C,  the  fact  that  the  strength  of  the 
C-C  t r  bond  is  nearly  twice  that  of  the  Si-C  n  bond  makes  the  second  resonance 
structure  much  more  important  This  accounts  for  the  fact  that  the  HjSi  =  CH  group 
has  such  a  large  effect  on  reducing  the  C-H  BDE.  Breaking  the  C-H  bond  in 
f^S^CH-CHj  allows  it  delocalization  to  occur,  resulting  in  the  Si-C  n  bond  in  the 
reactant  being  transformed  into  an  allylic  radical  in  which  the  it  bond  is  largely  C=C 


in  character. 
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The  same  allylic  radical  is  formed  by  breaking  the  Si-H  bond  in 
F^C^H-SiHg,  and  the  expected  localization  of  the  double  bond  in  the  radical  is 
evidenced  by  the  absence  of  any  appreciable  effect  of  the  F^C  =  CH  substituent  on 
the  Si-H  BDE  and  by  the  bond  lengths  in  the  radical,  which  are  nearly  the  same  as 
those  in  E^C=CH-SiH^.  The  predicted  inability  of  an  adjacent  C=C  n  bond  to 
provide  appreciable  delocalization  for  an  SiE^  center  is  consistent  with  the 
experimental  finding  that  the  SiH  BDE  is  unaffected  by  a  phenyl  substituent4  and 
that  of  the  EPR  spectrum  of  Ph-SiHj  shows  the  radical  to  be  largely  centered  at 
silicon.16  The  EPR  spectrum  of  [I^C=CH-SiH^]  likewise  has  the  radical  center 
primarily  localized  at  silicon  center,  although  there  is  some  delocalization  into  the 
C-C  it  bond.17 

For  Y  =  Z  =  Si,  and  A=C,  the  second  structure  in  equation  2.8  also  makes 
the  dominant  contribution  because  of  the  greater  strength  of  Si-C  compared  to  Si-Si 
it  bonds.11  However,  the  difference  between  Si-C  and  Si-Si  tt  bond  strengths  is 
smaller  than  that  between  Si-C  and  C-C  n  bonds,  thus  accounting  for  the  fact  that 
HjSi  =  SiH  has  a  smaller  effect  on  reducing  the  C-H  BDE  than  does  HjSi  =  CH. 

The  dominance  of  the  F^C  =  SiH-Sif^  resonance  structure  predicts  little 
stabilization  of  an  Si  radical  center  by  a  P^C  =  SiH  substituent.  However,  this  group 
effects  a  significant  reduction  in  the  Si-H  BDE,  as  does  a  P^Si  =  SiH  substituent.  In 
order  to  separate  possible  a  from  jt  effects  of  vinylic  substituent  groups  on  the  A-H 
BDEs  and,  particularly,  on  the  Si-H  BDEs,  once  again  the  effect  of  rotating  the 
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substituent  group,  X,  out  of  conjugation  in  X-AHj  was  calculated  for  each  substituent 
group. 

Section  2-5-3  Calculation  of  Allylic  Resonance  Energies 

Resonance  energies  in  radicals  containing  unsaturated  groups,  X,  are 
sometimes  estimated  by  the  reduction  in  the  A-H  BDE  caused  by  the  presence  of  X 
in  X-AH,.18  However,  due  to  the  existence  of  possible  a  effects  on  the  A-H  BDEs, 
the  energy  required  to  relate  an  unsaturated  X  group  out  of  conjugation  provides  a 
more  accurate  assessment  of  the  n  resonance  energies.  The  allylic  resonance 
energies  (AREs)  that  we  computed  in  this  manner  are  listed  in  Table  2-8. 


Table  2-8  CISD/6-31G*  Derived  Allylic  Resonance  Energies  (AREs)  for  Allylic 
Radicals,  [X-AHJ. 


AE  (keal/mol) 


Substituent 

Group  (X) 

A  =  C 

A  =  Si 

I^C=CH 

15.3 

0.2 

HjS^CH 

31.6 

0.8 

H^SiH 

8.2 

0.3 

HjSi=SiH 

13.6 

5.6 

HjS^CH4 

7.4 

f^Si  =  SiH* 

5.0 

‘All  centers  are  constrained  to  be  planar. 
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Section  2-5-3- 1  Unsymmetrical  Allylic  Radicals 

As  expected,  because  of  the  greater  strength  of  it  bonds  to  carbon  than  to 
silicon,  the  AREs  for  each  substituent  group  are  larger  when  the  partial  ir  bond  in 
the  allylic  radical  that  is  broken  leaves  the  unpaired  electron  on  carbon  rather  than 
silicon.  For  example,  in  [F^SiCHCFy  radical,  the  results  in  Table  2-8  predict  that 
it  should  cost  about  32  kcal/mol  to  rotate  about  the  C-C  bond,  but  rotation  about 
the  C-Si  bond  should  be  almost  without  cost.  The  difference  in  the  energy  required 
to  rotate  each  of  the  two  terminal  groups  out  of  conjugation  is  almost  exactly  equal 
to  the  difference  between  the  energies  of  the  C-C  (65  kcal/mol/  and  C-Si  (35 
kcal/mol /  it  bonds.  Similarly,  the  difference  of  13.3  kcal/mol  in  the  energies 
required  to  rotate  the  SiHj  and  CF^  groups  out  of  conjugation  in  [HjSiSiHCHj]  is 
equal  to  the  difference  in  Si-C  (35  kcal/mol)  and  Si-Si  (22  kcal/mol/  ir  bond 
energies. 

Section  2-S-3-2  Symmetrical  Allylic  Radicals 

When  the  terminal  groups  of  the  allylic  radical  are  identical,  rotation  of  either 
group  out  of  conjugation  must,  by  symmetry,  be  equivalent  to  rotation  of  the  other. 
Hence,  unlike  the  case  with  unsymmetrical  allylic  radicals,  comparison  of  the  AREs 
for  symmetrical  allylic  radicals  must  be  between  different  radicals. 

The  allyl  analogs  -  [HjC-SiH-CHJ*,  [HjSi-CH-SiFy*,  and  [H^i-SiH-SiHJ* 
-  like  the  parent  allyl  radical,  [F^C-CH-CHj]  ,  each  have  two  equivalent  resonance 
structures.  Simple  Huckel  theory  predicts  that  the  ARE  in  each  allylic  radicals  is 
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proportional  to  the  strength  of  the  A-Z  ir  bond  present  in  each  of  the  two  resonance 
structures.  As  a  result,  the  anticipated  trend  for  AREs  is  expected  to  be  : 
[HjC-CH-CIl,]*  >  [HjC-SiH-CH,]*  « [H,Si-CH-SiH,]*  >  [HjSi-SiH-SiH,]*.  More 
quantitatively,  the  ratio  between  the  AREs  should  be  approximately  9:5:5:3  based  on 
the  relative  strengths  of  the  x  bonds  that  are  involved.11 

Inspection  of  Table  2-8  reveals  that,  with  one  glaring  exception,  this  is  the 
ratio  of  the  calculated  AREs.  The  one  exception  is  a  computed  ARE  of  0.8 
kcal/mol  for  [H,Si-CH-Sif^]  ,  as  opposed  to  8.2  kcal/mol  computed  for 
[HjC-SiH-CHj]  .  As  discussed  above,  one  would  naively  expect  both  to  have  the 
same  ARE. 

However,  inspection  of  the  geometries  of  the  allylic  and  rotated-allylic 
structures  in  Appendix  1  discloses  significant  pyramidal ization  of  terminal  silicon 
groups  in  [HjSi-CH-Sify]  .  If  the  silicon  centers  are  constrained  to  remain  planar, 
the  ARE  for  [HjSi-CH-SiHj]  is  calculated  to  be  7.4  kcal/mol,  which  follows  the 
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anticipated  trend,  and  is  very  close  to  the  ARE  of  8.2  kcal/mol  calculated  for 
[HQ-SiH-CHj]  .  Due  to  the  fact  that,  unlike  silicon  centered  radicals,  those 
centered  on  carbon  tend  to  be  planar,  the  terminal  atoms  in  the  latter  allylic  radical 
are  planar  even  without  being  constrained  to  be  so. 

Inspection  of  the  geometries  also  reveals  that  [I^Si-SiH-SiHg]  contains  three 
pyramidal  silicon  centers.  However,  unlike  the  case  in  [F^Si-CH-SiHj]  ,  the  ARE 
does  not  appear  to  be  affected  by  pyramidalization,  since  constraining  the  silicon 
centers  to  be  planar  actually  decreases  the  calculated  ARE  from  5.4  kcal/mol  to  5.0 
kcal/mol.  This  finding  is  the  opposite  of  the  usual  result  that  constraining  atomic 
centers  to  planarity  increases  the  strength  of  the  tt  bonds  that  they  form,  since 
planarity  allows  better  overlap  of  the  orbitals  involved  in  n  bonding. 

Constraining  the  silicon  centers  in  the  [HjSi-SiH-SiHJ  system  to  planarity 
decreases  the  calculated  ARE  slightly,  since,  when  forced  to  be  planar,  the  preferred 
geometry  at  the  three  silicon  centers  is  artificially  altered  in  the  conjugated  geometry; 
whereas,  just  a  single  silicon  center  is  affected  in  the  rotated  geometry.  The  energy 
required  to  planarize  the  three  silicon  centers  in  the  conjugated  geometry  apparently 
exceeds  slightly  the  energy  required  to  planarize  the  non-conjugated  SiHj  center  in 
the  rotated  geometry. 

Section  2-5-4  Effects  of  Vinyl  Substituents  at  Unconjugated  Geometries 

As  discussed  in  Section  2-3,  A-H  BDEs  can  be  influenced  by  both,  a  and  v 
effects.  The  previous  section  describes  the  w  effects  of  vinyl  substituents  on  reducing 
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H^ZH-AI^  +  AH<  -  R^ZH-AH,  +  AH,  (2.10) 

the  A-H  BDEs;  the  remaining  component  is  a  donation/withdrawal  by  the  vinyl 
substituents.  This  effect  can  be  computed  from  the  energy  of  equation  2.10,  with  the 
singly  occupied  orbital  in  X-AE^  held  orthogonal  to  the  *  bond  in  X.  The  energies, 
which  are  reported  in  Table  2-9,  are  equal  to  the  difference  between  the  BDEs  in 
Table  2-7  and  the  AREs  in  Table  2-8. 

Table  2-9  a  Effects  vinyl  derivatives  on  the  CI-SD/6-31G*  derived  BDEs  using 
the  reaction1  [X-Aty*  +  AH«  -f  X-AHg  +  »AHg 


Vinyl  Substituent  Group  (X) 

A=C 

c/5 

n 

< 

H2C=CH 

1.9 

-0.9 

I^Si  =  CH 

4.7 

-1.0 

H^C=SiH 

4.9 

3.0 

R,Si  =  SiH 

5.3 

2.8 

“The  radical  orbital  in  the  reactant  is  orthogonal  to  the  k 

bond. 

Inspection  of  the  data  in  Table  2-9  reveals  that,  as  discussed  for  a  CHg 
substituent  in  Section  2-3-3,  the  Si-H  BDEs  are  very  slightly  increased  when  a  carbon 
is  bound  to  the  silicon  center;  but  are  reduced  when  a  silicon  is  the  a  atom  of  the 
substituent  group.19  Not  surprisingly,  a  electronic  effects  in  altering  the  Si-H  BDE 
arc  primarily  determined  by  the  substituent  atom  a  to  the  Si-H  bond  being  broken 
and  are  nearly  independent  of  the  atom  fi  to  it. 
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The  C-H  BDEs,  as  expected,  are  reduced  by  all  substituents.  What  is  not 
clear  at  first,  however,  is  why  there  is  a  nearly  3  kcal/mol  difference  between 
P^C  =  CH  and  HjSi  =  CH  substituents  on  reducing  the  C-H  BDE.  If  a  electronic 
effects  are  primarily  determined  by  the  a  atom  of  the  substituent  group,  then  these 
two  substituents  would  be  predicted  to  have  similar  effects  on  the  C-H  BDEs. 

As  discussed  in  Section  2-3-3,  carbon  centered  radicals  are  stabilized  by 
hyperconjugation.  The  Si-C  a  bond  in  H^Si  =  CH  is  weaker  than  the  C-C  a  bond  in 
HjC  =  CH;  and,  hence,  the  former  substituent  might  hyperconjugate  better  than  the 
latter.  However,  this  explanation  cannot  be  verified  directly  with  these  vinyl 
substituent  since  rotation  of  the  Si-C  a  bond  out  of  conjugation  with  the  radical 
orbital  would  result  in  alignment  of  the  n  bond  with  the  radical  orbital.  Comparison 
of  the  calculated  BDE  for  SiHgCHjCHj-H  with  that  for  Cf^-CH^-H  can  elucidate 
if  a  carbon  centered  radical  can  more  effectively  hyperconjugate  with  a  C-Si  a  bond 
than  with  C-H.  Rotation  of  the  Si-C  a  bond  out  of  conjugation  with  the  radical 
center  in  HgSi-CHj-CHj  provides  another  test  of  the  effect  of  hyperconjugation  with 
an  Si-C  a  bond  on  stabilization  of  a  carbon  centered  radical. 

The  calculated  CI-SD  BDE  of  HgSiCHjCHj-H  is  2.8  kcal/mol  less  than  that 
of  CHjCHj-H  when  the  Si-C  a  bond  is  aligned  with  the  orbital  containing  the 
unpaired  electron  in  the  radical.  Rotation  about  the  C-C  bond  by  90°  in  the 
HgSiCHjCHs  radical  requires  2.2  kcal/mol  showing  that  the  lowering  of  the  BDE  is 
largely  due  to  Si-C  a  bond  hyperconjugation.  Therefore,  the  a  effect  calculated  for 
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the  F^S^CH  substituent  on  lowering  the  C-H  BDE  probably  also  actually  includes 
a  contribution  from  delocalization  of  the  carbon  radical  through  hyperconjugation. 

Hyperconjugation  with  a  C-Si  bond  might  also  be  expected  to  be  important 
in  stabilizing  [F^C-SiH-CHj]  when  the  n  bond  of  the  vinyl  group  is  twisted  out  of 
conjugation  with  the  singly  occupied  orbital  on  carbon.  However,  since  Si-C  n  bonds 
are  weaker  than  C-C  n  bonds,  hyperconjugation  should  be  less  important  in  twisted 
[HjC-SiH-CHg]  .  In  order  to  test  this  prediction,  we  performed  calculations  on  the 
BDE  of  CHgSiHgCH^-H  and  the  barrier  to  rotating  the  Si-C  a  bond  out  of 
conjugation  in  the  resulting  radical. 

The  calculated  CI-SD  BDE  of  CE^SiF^CHj-H  is  1.0  kcal/mol  weaker  than 
that  of  HgSiCH^H,  and  the  barrier  to  twisting  the  Si-C  bond  out  of  conjugation  in 
CHjSiHjCHj  is  only  0.1  kcal/mol.  Clearly,  stabilization  of  this  radical  by 
hyperconjugation  with  the  Si-C  a  bond  is  less  important  than  in  SiF^Cf^CHj. 

This  finding  implies  that  the  similar  C-H  BDEs  for  non-conjugated  HjSi  =CH 
and  I^C  =  SiH  substituents  is  due  to  the  fact  that,  although  the  a  bond  in  the  former 
substituent  hyperconjugates  better,  the  latter  has  the  more  electropositive  silicon 
atom  in  the  a  position.  Further  evidence  for  the  relative  unimportance  of  a  bond 
hyperconjugation  when  silicon  is  in  the  a  position  of  the  substituent  group  comes 
from  the  finding  in  Table  2-9  that  a  twisted  F^S^SiH  substituent  lowers  the  C-H 
BDE  only  0.4  kcal/mol  more  than  a  twisted  H2C  =  SiH  substituent;  whereas  with 
carbon  in  the  a  position,  twisted  F^Si  =  CH  lowers  the  C-H  BDE  by  2.8  kcal/mol 
more  than  twisted  11^= CH. 
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Section  2-6  Conclusions 

As  noted  in  section  2-3-4,  we  found  that  a  substituent  affects  A-H  BDEs  in 
two  ways.  First,  the  well  understood  effect  of  rr  delocalization  in  X-AF^  wealems 
A-H  bonds  by  stabilizing  the  radicals.  Second,  a  electron  donors  weaken  A-H  bonds, 
whereas  a  electron  withdrawing  substituents  strengthen  them. 

Both  of  these  effects  play  an  important  role  in  determining  the  strength  of 
C-H  bonds,  rr  stabilization  of  carbon  centered  radicals  usually  can  overcome  the 
effects  of  a  electron  withdrawal,  except  when  a  withdrawl  is  severe,  as  in  FSCH. 

Conversly,  the  Si-H  bond  is  most  strongly  influenced  by  a  electronic  effects. 
Of  all  the  substituents  examined,  only  two  --  BHj,  I^Si  =  SiH  -  have  been  shown  to 
provide  a  significant  amount  of  rr  delocalization. 
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bond,  but  the  resultant  it  bond  that  is  formed  is  weaker  than  the  it  bond  formed 
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energies  of  the  it  bonds  formed  to  SiHj  and  CHj  substituents  is  almost  exactly  equal 
to  the  difference  in  the  energies  of  the  SiH  and  CH  bonds  being  broken.  For  a 
carbon  centered  radical,  CHg  should  be  slightly  more  stabilizing  than  SiHj. 
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13.  The  increase  in  C-H  BDE  from  100.3  kcal/mol  for  HgCCHj,  to  106  kcal/mol  for 
HjCCHj  and  to  132  kcal/mol  for  HCCH  demonstrates  this  fact  very  well.  All  values 
are  from:  CRC  Handbook  of  Chemistry  and  Physics,  67th  ed.;  Weast,  R.  C.,  Ed.;  CRC 
Press:  Boca  Raton,  FL>  1986;  p  F-178. 

14.  For  leading  references,  see  Pasto  (1). 

15.  Borden,  W.  T.  "Modern  Molecular  Orbital  Theory  for  Organic  Chemists," 
Prentice-Hall,  Inc.,  1975. 

16.  Sakurai,  H.;  Umino,  H.;  Sugiyama,  H.  J.  Am.  Chem.  Soc.,  1980,  102,  6837. 

17.  Jackson,  R.  A.;  Zarkadis,  A.  K.  Tetrahedron  Letters,  1988,  29,  3493. 

18.  See,  for  example,  Doering,  W.  v.  E.;  Roth,  W.  R.;  Breuchkmann,  R.;  Figge,  L.; 
Lennartz,  H.-W.;  Printzbach,  H.  Chem.  Ben,  1988, 121 ,  1. 

19.  The  values  in  Table  2-7  are  calculated  at  the  CI-SD/6-31G*  level  of  theory. 
Since  the  BDE  values  given  in  Section  2-3-3  were  calculated  using  the 
MP4SDTQ/6-31G*  level  of  theory,  these  values  were  recomputed  at  the 
CI-SD/6-31G*  level  for  direct  comparison.  The  CI-SD  values  for  the  difference  in 
C-H  BDEs  are  4.3  kcal/mol  for  ethane,  5.6  kcal/mol  for  silylmethane.  The 
difference  in  the  Si-H  BDEs  are  -0.1  kcal/mol  for  methylsilane,  and  3.5  kcal/mol  for 
silylsilane.  The  CI-SD  values  result  in  approximately  1.0  kcal/mol  additional 
stabilization  over  MP4SDTQ  values  for  the  difference  in  BDEs  of  the  C-H  bond  and 
0.6  kcal/mol  for  the  Si-H  bond.  The  difference  in  energy  is  a  direct  result  of  the 
wave  function  of  the  more  substituted  radical  containing  a  greater  amount  of  spin 
contamination  from  the  quartet  state.  M011er-Plesset  perturbation  theory  then 
spurriously  overestimates  the  BDE  of  the  more  highly  substituted  radical,  resulting 
in  a  difference  of  BDEs  that  is  too  small. 


Chapter  3  Pyramidalization  of  Phosphorus  in  Borylphosphines 


Section  3-1  Chemical  Background 
Section  3-1-1  Phosphorus  versus  Nitrogen 

Trivalent  phosphorous  compounds  differ  from  their  nitrogen  analogs,  not  only 
by  having  longer  bond  lengths  but  also  by  being  highly  pyramidal1  and  having  much 
larger  barriers  to  inversion.2  Indeed,  unlike  their  nitrogen  analogs,  optically  active 
phosphines  can  be  prepared  and  racemize  only  at  elevated  temperatures.3 

Due  to  the  highly  pyramidal  geometry  and  large  energies  required  for 
planarization  of  phosphines,  the  lone  pair  on  phosphorous  is  much  less  available  for 
7r  bonding  than  the  lone  pair  of  amines.4,5  For  example,  phosphines  are  not  planar, 
even  when  attached  to  three  formyl  groups;  and  the  barrier  to  rotation  about  the 
phosphorus-carbonyl  bond  is  small;  whereas  amides  have  planar  geometries  at 
nitrogen  and  considerably  larger  barriers  to  rotation  about  the  nitrogen-carbonyl 
bond. 

Section  3-1-2  Competition  for  Electron  Density 

The  rotational  barrier  about  the  nitrogen  carbonyl  bond  in  an  amide  decreases 
when  multiple  carbonyls  are  attached  to  a  central  nitrogen.  When  only  one  formyl 
group  is  attached,  the  rotational  barrier  about  the  C-N  bond  is  17.8  kcal/mol. 
However,  the  barrier  goes  down  to  12.9  kcal/mol  with  the  addition  of  a  second 
formyl  group,  and  a  third  formyl  reduces  the  rotational  barrier  all  the  way  to  7.5 
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kcal/mol.6  Apparently,  the  formyl  groups  compete  for  the  electron  density  in  the 
nitrogen  lone  pair  orbital,  which  is  why  n  bond  strength  and,  hence,  the  rotational 
barrier  decreases  with  increasing  number  of  formyl  groups. 

The  same  type  of  competition  has  been  computationally  predicted  to  exist  in 
aminoboranes,  borylamines,  hydroxyboranes  and  borylhydroxides.7  Both  of  these 
types  of  compounds  have  lone  pairs  of  electrons  on  atoms  that  compete  to  donate 
electron  density  into  an  empty  boron  n  orbitals.  In  both  cases,  the  greatest  barrier 
to  rotation  occurs  when  there  is  just  one  substituent  with  a  lone  pair  adjacent  to  a 
single  boron  with  its  empty  orbital.  When  there  is  an  excess  of  either  boron 
substituents  attached  to  a  center  with  a  lone  pair  or  substituents  with  lone  pairs 
attached  to  a  boron  center,  there  is  reduction  in  barrier  to  rotation.8 

Competition  effects  are  also  found  in  the  rotational  barriers  about  C-C  bond 
along  the  series  (CH^CHg.,,,  where  n  =  l  is  ethene,  n=2  is  allyl,  and  n=3  is 
trimethylenemethane.  The  strongest  ir  bond  is  that  in  ethene,  where  rotation,  which 
breaks  a  full  tt  bond,  requires  65  kcal/mol.9  In  allyl,  the  barrier  to  rotation  of  one 
methylene  is  reduced  to  15  kcal/mol.10  Calculations  reproduce  these  two  barriers 
rather  well.11  In  triplet  trimethylenemethane,  the  barrier  height  has  not  been 
measured;  but  it  has  been  calculated  to  be  10  kcal/mol  at  comparable  levels  of 
theory.12 

Competition  effects  on  the  barriers  to  C-C  rotation  in  the  series  (CH^CHj.,, 
are  predicted  qualitatively  even  by  simple  Huckel  molecular  orbital  theory.  The 
HMO  energies  of  the  bonding  MOs  are  respectively  a  +  /3  for  n=l,  a  +  1.414  0 
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for  n=2,  and  a  +  1.732  0  for  n=3.ls  Hence,  the  energy  increase  on  breaking  one 
n  bond  in  each  of  these  hydrocarbons  is  2  /3,  0.828  (3,  and  0.636  /?  respectively.14 

Given  the  apparently  ubiquitous  nature  of  competitive  substituent  effects,  the 
discovery  of  a  cooperative  substituent  effect  would  be  noteworthy.  As  discussed  in 
Section  3-3,  we  discovered  such  an  effect  in  HP(BI^)2  and  performed  calculations 
to  establish  its  origin. 

Section  3-1-3  Experimental  Examples  of  Electron  Delocalization  from  Phosphorus 

Despite  the  fact  that  phosphorus  has  a  high  barrier  to  planarity,  there  is 
experimental  evidence  that  in  boryl  phosphines,  the  empty  boron  n  orbital  is  effective 
in  delocalizing  the  phosphorous  lone  pair.  The  x-ray  structure  of  a  monoboryl 
phosphine,  (diphenylphosphino)dimesitylborane,  shows  a  pyramidalization  angle,  <p 
(the  angle  between  the  extension  of  the  P-B  bond  and  the  plane  containing 
phosphorus  and  the  other  two  groups  bound  to  it)  of  43° 15  compared  to 
pyramidalization  angles  of  about  68oie  in  the  unconjugated  phosphine,  PPhj.  In  a 
similarly  substituted  derivative  of  diboryl  phosphine,  the  boron-phosphorous  bond 
lengths  are  shorter17  and  the  phosphorous  is  now  planar  (0  =  0°  ).18 

There  is  also  evidence  for  strong  P-B  n  bonding  in  the  crystal  structures  of 
cyclic  phosphorous  boron  compounds.  The  six-membered  ring  in  (mesityl  B-P 
cyclohexyl)^  is  planar  and  exhibits  equal  P-B  bond  lengths  of  1.84A.19  The  four 
membered  ring  in  (2,2,6,6-tetramethylpiperidio  B-P  2,4,6-tris-t-butylbenzene)2  is  not 
planar;20  but  this  is  not  surprising,  since  planarity  would  yield  a  conjugated  cyclic 
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4-electron  rr  system.  Huckel  MO  theory  predicts  that  this  it  system  would  be 
antiaromatic;  hence,  delocalization  of  the  phosphorus  lone  pairs  would  be  less 
stabilizing  than  in  the  molecule  containing  a  six-electron  cyclic  it  system. 

In  order  to  examine  computationally  the  effect  of  boron  substituents  on  the 
barrier  to  inversion  of  phosphorous,  we  performed  calculations  on  (BH^PHg_n, 
n=l-3,  and  on  the  cyclic  (HBPH)^  compounds  with  n  =  2  and  n  =  3. 

Section  3-2  Theoretical  Methodology 

Calculations  were  carried  out  using  Gaussian  8221  and  8< 92  with  the  6-31G* 
basis  set.  All  geometries  were  first  optimized  using  a  RHF  wave  function.  Electron 
correlation  was  added  at  these  geometries  using  fourth  order  Miller  Plesset 
perturbation  theory,  MP4SDTQ.  Some  geometries  were  reoptimized  using  MP2,  to 
verify  that  inclusion  of  electron  correlation  during  optimization  did  significantly  alter 
the  geometry  or  the  MP4SDTQ  energy  of  the  RHF  optimized  molecule.  Zero  point 
energy  corrections  from  RHF  vibrational  analyses  were  included. 

The  four  geometries  optimized  for  H^PBF^  are  the  same  as  those  chosen  by 
Groperf3  for  his  RHF  calculations.  The  geometries  are  pictured  in  Figure  3-1. 
They  are:  1,  the  equilibrium  Q  geometry  with  the  lone  pair  and  empty  orbital  in  the 
symmetry  plane;  2,  the  completely  planar,  Qv,  structure;  3,  the  rotated  Q  structure 
with  the  boron  hydrogens  in  the  symmetry  plane  and  the  H-P-H  bond  angle  bisected 
by  the  BP^  group;  and  4,  the  rotated  Qv  structure  with  a  planar  phosphorus. 
Reoptimization  of  the  Q  equilibrium  (1)  and  the  completely  planar  Qv  structure  (2) 
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at  the  MP2  level  resulted  in  a  decrease  in  the  central  bond  length  of  l-2%.2'*  The 
change  in  the  MP4SDTQ  energies  relative  to  each  other  amounted  to  less  than  0.5 
k cal /mol 25  These  differences  were  not  considered  large  enough  to  warrant  using 
electron  correlation  during  optimization  of  any  of  the  other  structures. 

For  HPCBHjX,  we  optimized  the  six  structures  in  Figure  3-2.  They  are:  5,  the 
equilibrium  Q  structure  with  both  empty  boron  orbitals  aligned  with  the  phosphorous 
lone  pair;  6,  the  completely  planar  Qv  structure;  7,  the  transition  state  for  rotation 
of  one  BHj  out  of  conjugation  with  the  phosphorus  lone  pair;  8,  the  Q  structure  with 
one  BF^  rotated  out  of  conjugation  with  a  planar  phosphorus;  9,  the  dirotated  Q 
transition  state  structure;  and  10,  the  dirotated  C,h  structure  with  a  planar  geometry 
at  phosphorous.  The  equilibrium  (5)  and  the  completely  planar  (6)  structures  were 
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Figure  3-2  geometries. 


again  reoptimized  with  inclusion  of  electron  correlation  at  the  MP2  level;  and,  once 
again,  little  change  was  found  in  the  optimized  geometries  or  MP4SDTQ  energies.26 

Vibrational  analyses  were  performed  at  all  stationary  points.  The 
monorotated  geometries,  3  and  7,  were  found  to  have  exactly  one  imaginary 
frequency,  corresponding  to  rotation  of  the  twisted  boron  back  into  conjugation,  and 
each  was  thus  demonstrated  to  be  true  transition  states.  In  contrast,  the  dirotated 
Q  structure  (9)  was  found  to  have  two  imaginary  frequencies.  Hence  it  is  not  a  true 
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saddle  point  but  is  only  a  transition  state  when  constrained  by  symmetry.  One 
vibrational  mode  with  an  imaginary  frequency  does,  in  fact,  correspond  to  disrotary 
motion  of  the  two  groups,  which  preserves  the  symmetry  plane;  but  the  other 
corresponds  to  conrotary  motion,  which  destroys  the  symmetry  plane. 

Five  geometries  of  triborylphosphines  were  studied.  These  geometries,  which 
are  depicted  in  Figure  3-3,  were:  11,  the  equilibrium,  Dsh,  planar  structure  in  which 
the  empty  p  orbitals  on  each  boron  is  aligned  for  overlap  with  the  phosphorous  lone 
pair;  12,  the  non-planar,  monorotated,  Q  structure;  13,  the  planar,  monorotated,  Qv 
structure;  14,  the  planar,  dirotated,  Qv  structure;  15,  the  planar,  trirotated,  D3h 
structure  in  which  the  empty  boron  orbitals  cannot  conjugate  with  the  phosphorous 
lone  pair.  No  attempt  was  made  to  optimize  any  of  these  geometries  with  the 
inclusion  of  electron  correlation  at  the  MP2  level. 

Three  conformations  of  the  four  membered  ring,  P^Ff,,  were  examined. 
These  were  the  completely  planar  structure,  16,  a  puckered  Qv  structure,  17,  with  the 
hydrogens  at  phosphorus  pyramidalized  cis  to  each  other,  and  a  Qh  structure,  18, 
with  the  phosphorous  hydrogens  pyramidalized  trans  to  each  other  and  the  rest  of  the 
molecule  coplanar.  For  the  6-membered  ring  the  conformations  examined  were,  the 
planar  D3h  benzene-like  geometry,  and  the  non-planar,  Qh  geometry,  20,  which 
resembles  chair  cyclohexane  with  the  phosphorus  lone  pairs  taking  the  place  of  three 
axial  hydrogens.  Since  these  cyclic  systems  were  rather  large  and  required  a 
significant  amount  of  cpu  time  to  optimize  at  even  the  RHF  level,  we  did  not  try 
reoptimizing  any  of  them  at  the  MP2  level. 
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Figure  3-3  P(BHj)^  geometries 


Figure  3-4  P2Bl,H,  and  PjB^  geometries. 
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Section  3-3  Results 

Section  3-3-1  Monoboryl  Phosphine 


Table  3-1  Relative  Energies  (kcal/mol)  for  P^PBF^  Calculated  with  the  6-31G* 


s  Set. 

Geometry 

Svmmetrv 

RHF 

MP4SDTO 

Equilibrium  (1) 

Q 

0 

0 

Planar  (2) 

C2  V 

8.4 

4.2 

Rotated  (3) 

Q 

7.9 

9.8 

Planar  Rotated  (4) 

Qv 

41.9 

43.3 

A  phosphorus  with  a  single  BHj  substituent  was  calculated  to  have  a 
pyramidal  equilibrium  geometry.  This  finding  is  in  agreement  with  earlier 
calculations?3 .and  the  experimental  jesultSL  on  a  highly  substituted.derivativa.15 .  Tbe  . 
calculated  pyramidalization  angle,  <p,  at  phosphorus  of  71  °  in  HjPBP^  is  less  than 
that  of  86°  in  PHg;27  and  the  MP4SDTQ  barrier  to  planarization  of  4.2  kcal/mol 
in  I^PBHj  is  substantially  reduced  from  the  35.8  kcal/mol  calculated  for  PH3  .  At 
the  MP4SDTQ  level,  it  requires  about  half  the  energy  to  planarize  the  phosphorous 
as  it  does  to  rotate  the  boron  out  of  conjugation.  This  last  prediction  is  the  opposite 
of  that  made  by  earlier  RHF  calculations  by  Gropen.23  He  found  that  rotation  of  the 
boron  is  easier  than  inversion  of  the  phosphorous  --  6.6  kcal/mol  versus  8.1  kcal/mol. 
Although,  his  calculations  did  not  include  d  polarization  functions  on  the  boron,  a 
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larger  barrier  to  inversion  than  to  rotation  is  also  found  by  our  RHF  calculations 
with  the  polarized  6-31G*  basis  set. 

The  inclusion  of  electron  correlation  favors  n  electron  delocalization.28  This 
decreases  the  barrier  to  planarization  of  the  phosphorus  and  increases  the  barrier  to 
rotation  of  the  boron  out  of  conjugation.  Planarization  of  the  phosphorus  with  the 
empty  boron  orbital  rotated  out  of  conjugation  requires  34.0  kcal/mol  at  the  RHF 
level  and  33.5  kcal/mol  at  the  MP4SDTQ  level,  which  are  close  to  the  values  of  37.8 
kcal/mol  and  35.8  kcal/mol  for  PHg  at  the  same  levels  of  theory.  Basis  set 
differences  apparently  account  for  the  difference  between  our  RHF  value  for 
planarization  of  HjPBRj  and  Gropen’s  value  of  42.0  kcal/mol. 

The  effect  of  a  single  boron  group  conjugated  with  the  phosphorous  lone  pair 
is  not  enough  to  planarize  the  phosphine,  but  it  does  make  planarization  much 
easier.  Indeed,  when  conjugated  with  the  BHg  group,  the  phosphorus  resembles  an 
unconjugated  amine  in  its  barrier  to  inversion.  The  stabilization  through  tt 
delocalization  increases  as  phosphorus  becomes  more  planar,  since  at  a  planar 
geometry  the  tt  overlap  between  the  phosphorus  lone  pair  and  the  empty  orbital  on 
boron  is  maximized.  This  accounts  for  the  high  barrier  to  rigid  rotation  of  the  BHj 
group  when  phosphorus  is  held  planar. 

Section  3-3-2  Diboryl  Phosphine 

With  two  BH,  groups,  the  potential  energy  surface  becomes  more 
complicated.  In  addition,  we  find  that  predicting  the  equilibrium  geometry'  accurately 
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becomes  difficult  because  the  potential  surface  for  planarization  at  phosphorus  is 
remarkably  flat.  The  MP2  optimized  pyramidalization  angle  at  phosphorous  is 
40°, 29  but  the  barrier  to  planarity  is  only  0.3  kcal/mol  at  the  MP4SDTQ  level. 
Thus,  the  hydrogen  on  the  phosphorus  can  be  swung  through  more  than  80°  of  arc30 
with  hardly  any  change  in  energy.  Although  the  minimum  on  the  potential  surface 
occurs  at  a  significantly  pyramidalized  geometry,  the  experimental  finding  that  a 
highly  substituted  derivative  of  HP^H^  is  planar31  is  not  at  all  surprising. 
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Table  3-2  Relative  Energies  (kcal/mol)  for  I  IP(BIT2)t2  Calculated  Using  the 


■31G*  Basis  Set 

Geometry 

Svmmetrv 

RKF 

MP4SDTO 

Equilibrium  (5) 

Q 

0 

0 

Planar  (6) 

Qv 

0.4 

0.3 

Monorotated  (7) 

Q 

13.1 

16.6 

Planar  Monorotated  (8) 

Q 

19.9 

20.2 

Dirotated  (9) 

Q 

20.4 

24.8 

Planar  Dirotated  (10) 

Qv 

53.5 

59.5 

One  very  surprising  result  of  the  calculations  on  the  HP^H^  surface  is  that 
rotation  of  the  first  BP^  group  out  of  conjugation  requires  16.6  kcal/mol  whereas 
rotation  of  the  BHj  group  in  P^PBH,  required  about  half  that,  9.8  kcal/mol. 
Rotation  of  the  second  BHj  group  out  of  conjugation  in  HP^H^  requiring  8.2 
kcal/mol  is  consistent  with  the  value  in  HjPBHj.  Since  the  first  BHj  group  in 
HP(BH2)2  requires  more  energy  to  rotate  out  of  conjugation  than  the  second,  it 
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would  appear  that  the  empty  orbitals  on  the  two  borons  do  not  compete  for  electron 
density  from  the  phosphorous  lone  pair,  but  instead,  somehow  act  cooperatively. 

Never  the  less,  the  empty  orbitals  on  each  boron  can  be  shown  xo  compete  for 
phosphorus  electron  density  when  the  phosphorus  is  constrained  to  remain  planar. 
With  this  constraint,  rotation  of  the  first  BHg  group  out  of  conjugation  requires  19.9 
kcal/mol;  and  rotation  of  the  second  requires  an  additional  39.3  kcal/mol,  which  is 
very  close  to  the  value  of  39.1  kcal/mol  for  rotation  in  planar  HgPBFLj.  In  planar 
HP(BH2)2,  competition  between  the  two  BF^  groups  apparently  makes  the  energy 
required  for  rotation  of  the  first  BHg  out  of  conjugation  half  of  that  for  the  second. 
Once  the  first  BHj  group  has  been  rotated  out  of  conjugation,  there  is  virtually  no 
effect  of  the  first  BFL,  group  on  rotation  of  the  second  BHg  group  out  of  conjugation. 

Based  on  the  calculated  barriers  to  BF^  group  rotation,  two  BHj  groups,  both 
conjugated  with  the  lone  pair  on  a  planar  phosphorus,  compete  for  electron  density; 
but,  when  the  phosphorus  is  allowed  to  pyrainidalize,  the  two  BF^  groups  appear  to 
cooperate  with  one  another  in  delocalizing  the  phosphorus  lone  pair.  Thus, 
depending  on  whether  or  not  phosphorus  is  constrained  to  remain  planar,  our 
calculations  find  that  the  two  borons  can  cooperate  with  each  other.  What  happens 
with  the  addition  of  a  third  BF^  group? 

Section  3-3-3  Triboryl  Phosphine 

According  to  our  calculations,  the  equilibrium  structure  about  the  trivalent 
phosphorous  in  P(BF^)j  is  planar.  This  result  was  confirmed  by  a  frequency  analysis 
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which  gave  no  imaginary  frequencies,  so  the  planar  geometry  is  a  true  minimum. 
The  potential  surface  for  pyramidalization  of  monorotated  P(Bf^)^  is  as  flat  as  that 
of  fully  conjugated  HP(BH^  surface.  The  equilibrium  geometry  occurs  at  <p  = 
25.5°,  but  the  barrier  to  inversion  is  only  0.3  kcal/mol. 


Table  3-3  Relative  Energies 

(kcal/mol)  for  P(BI^^  Using  the  6-3 1G 

Geometry 

Svmmetrv 

RHF 

MP4SDTO 

Equilibrium  (11) 

^3h 

0 

0 

Monorotated  (12) 

Q 

14.4 

12.8 

Planar  Monorotated  (13) 

Qv 

14.4 

13.0 

Planar  Dirotated  (14) 

Qv 

34.0 

32.8 

Planar  Trirotated  (15) 

^3h 

68.1 

72.7 

Rotation  of  a  single  BJ-^  group  out  of  conjugation  in  P(Bf^)^,  requires  12.8 
kcal/mol,  which  is  less  than  the  16.6  kcal/mol  in  HPtBH^.  Thus,  it  appears  that 
the  three  BHa  groups  in  PtBH^  are  competing  for  electron  density,  instead  of 
cooperating.  The  competition  is  also  seen  in  P^H^  when  rotating  BE^  groups  out 
of  conjugation  while  maintaining  a  planar  phosphorus  center.  The  energy  required 
for  rotation  of  the  first  is  13.0  kcal/mol;  the  barrier,  for  rotation  of  the  second  BHj 
group  is  19.8  kcal/mol;  and  rotation  of  the  third  BHa  group  requires  39.9  kcal/mol. 
Once  again,  the  latter  two  values  are  almost  identical  to  those  calculated  for  rotation 
about  the  P-B  bond,  while  maintaining  a  planar  phosphorus  in,  respectively, 
HP(BHa)a  and  HjPBHj. 
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Section  3-4  Competition  versus  Cooperation 

Based  on  competitive  substituent  effects,  the  anticipated  trend  for  rotation  of 
one  BHj  group  out  of  conjugation  in  Hj.JPfBP^  would  be  progressively  lower 
barriers  to  rotation  as  the  number  of  BHj  groups  bound  to  a  central  phosphorus 
increases.  However,  as  shown  by  the  data  summarized  in  Table  3-4,  RjPBHj  has  a 
lower  barrier  than  either  HP^H^  or  P(BE^)j.  This  finding  suggests  the  existence 
of  a  second  substituent  effect  that  is  cooperative,  rather  than  competitive.  This 
second  effect  appears  to  vanish  when  the  phosphorus  is  constrained  to  be  planar, 
since  the  predicted  barriers  to  rotation  are  restored  to  the  order  anticipated  on  the 
basis  of  pure  competition  between  substituents. 

Table  3-4.  MP4SDTQ  Barrier  to  Rotation  of  One  BH,  Group  Out  of 
Conjugation  in  H3_JP(Bf^)^,  n  =  l-3,  With  and  Without  the  Constraint  of  a 
Planar  Phosphorus. 


AE^t  (kcal/mol) 


Compound 

4>conjS 

<pj 

Equilibrium 

Structures 

Planar 

Phosphorus 

f^PBHj 

71.0° 

82.4' 

9.8 

39.1 

HPfBEU 

45.1' 

71.9* 

16.6 

19.9 

PCBH,)* 

0.0' 

25.5' 

12.8 

13.0 

“Pyramidalization  angle  at  phosphorus  in  conjugated  structure. 
bPyramidalization  angle  at  phosphorus  in  rotated  structure. 
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Figure  3-5  Separation  of  the  Total  Energy  into  Planarization  Energy  and 
Conjugation  Energy. 


As  summarized  in  Table  3-4,  there  are  changes  in  the  geometry  at  the 
phosphorus  center  that  occur  as  more  BH,  groups  are  added,  namely,  phosphorus 
becomes  more  planar  (i.e.  <p  gets  smaller).  There  is  also  a  change  in  geometry  at 
phosphorus  that  occurs  on  rotation  of  a  group  out  of  conjugation.  In  order  to 
separate  the  changes  in  energy  that  accompany  the  changes  in  geometry  at 
phosphorus  from  the  pure  conjugation  energy  of  each  BHj  group,  we  examined  the 
potential  energy  surfaces  of  HjPBHj  and  HP^H^  in  greater  detail  using  the 
energies  of  the  reactions  depected  in  Figure  3-5.  Since  two  sets  of  calculations  were 
performed  at  each  of  several  values  of  the  pyramidalization  angle,  <p ,  these 
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calculations  were,  for  the  sake  of  economy,  carried  out  at  the  RHF  rather  than 
MP4SDTQ  level  of  theory. 

Section  3-4-1  The  H2PBH2  Potential  Surface 

The  equilibrium  structure  of  H^PBH^  contains  a  pyramidal  phosphorus  and 
a  rather  weak  n  bond  between  the  phosphorus  and  boron.  Since  the  energy  required 
to  rotate  the  BHj  group  increases  substantially  when  the  phosphorus  is  constrained 
to  be  planar,  the  7r  bond  strength  apparently  increases  as  the  phosphorus  becomes 
more  planar.  This  is  depicted  graphically  in  Figure  3-6. 

The  RHF  conjugation  energy  --  the  negative  of  the  energy  required  to  rotate 
the  BHj  group  out  of  conjugation  with  the  phosphorus  lone  pair  --  is  plotted  as  a 
function  of  <p.  Also  plotted  is  the  RHF  planarization  energy  --  the  energy  required 
to  flatten  the  phosphorus  with  the  BP^  group  twisted  out  of  conjugation,  relative  to 
the  RHF  energy  of  the  monorotated  equilibrium  geometry,  which  has  0  =  82.4 0 . 
The  sum  of  these  two  energies  is  the  RHF  total  energy  -  the  RHF  energy  of  the 
conjugated  molecule,  relative  to  the  RHF  energy  of  the  optimal,  monorotated 
structure.  The  total  energy  curve  in  Figure  3-6  is  that  which  is  traversed  when  the 
phosphorus  center  inverts  in  conjugated  HjPBHj. 

Since  the  conjugation  energy  and  planarization  energy  both  have  a  functional 
dependance  on  0  that  is  similar  in  magnitude  but  opposite  in  sign,  addition  of  these 
two  energies  gives  a  total  energy  curve  that  is  rather  flat.32  Therefore,  although  a 
conjugated  BF^  group  reduces  the  equilibrium  value  of  0  only  slightly,  from  82°  to 
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Figure  3-6  One  dimensional  P^PBHj  potential  energy  surface  in  <p. 


71  *,  conjugation  reduces  the  barrier  to  planarity  substantially,  from  34.0  kcal/mol  to 
8.4  kcal/mol  at  RHF  level  and  from  33.5  kcal/mol  to  4.2  kcal/mol  at  MP4  level.33 


Section  3-4-2  The  HP(BH2)2  Potential  Surface 

As  discussed  above,  in  HP(BH,)^  the  barrier  to  rotation  of  one  BHj  out  of 
conjugation  appears  to  indicate  that  the  two  BH,  substituents  are  cooperating,  rather 
than  competing.  An  important  clue  to  the  origin  of  this  apparent  cooperative  effect 
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is  the  shape  of  the  curve  in  Figure  3-7  for  the  planarization  energy.  This  curve 
resembles  more  closely  the  curve  in  Figure  3-6  for  the  total  energy  of  HjPBHj 
instead  of  the  planarization  energy  of  HjPBHj. 

Since  the  potential  curve  for  planarizing  monorotated  HP(BHj),  is  relatively 
flat,  <p  can  be  decreased,  in  order  to  form  strong  P-B  >r  bonds  in  the  fully  conjugated 
structure,  at  a  fraction  of  the  cost  required  in  rotated  HjPBHj.  Because  the  curve 
for  the  planarization  energy  in  HPCBHj),  is  quite  flat,  the  total  energy  curve  is  flatter 


still.  Indeed,  at  the  RHF  level,  only  0.5  kcal/mol  is  required  to  go  from  the 
equilibrium  structure,  at  0  =  45%  to  the  planar  structure. 

Although  the  softer  potential  surface  for  planarization  of  monorotated 
HP(BH2)2  results  in  an  apparent  cooperative  effect  between  the  two  BF^ 
substituents;  comparison  of  the  curves  for  the  conjugation  energies  in  Figure  3-6  and 
Figure  3-7  reveals  the  expected  competition  in  HP(BFlj)^  between  the  two  BHj 
groups.  Hence,  at  most  values  of  <t> ,  the  magnitude  of  the  conjugation  energy  (the 
energy  required  to  twist  a  BFLj  group  out  of  conjugation  at  fixed  <p)  is  larger  in 
HjPBF^  than  HP^FLj^.  For  example,  as  noted  earlier,  when  the  phosphorus  is 
constrained  to  be  planar,  competitio  l  is  seen  (Table  3-4);  and  the  barrier  to  rotate 
a  BF^  group  out  of  conjugation  in  Hg.^BH^  decreases  as  n  increases. 

However,  at  large  values  of  <p,  careful  comparison  of  Figures  5  and  6  shows 
that  two  BHj  groups  do  cooperate,  since  the  conjugation  energy  in  HP(BF^  is 
slightly  larger  than  that  in  FLjFBF^  for  <p  >  60°.  At  these  highly  pyramidalized 
geometries,  the  phosphorus  lone  pair  has  a  large  amount  of  s  character  and  is  not 
close  to  being  parallel  to  the  empty  orbitals  on  boron.  Indeed,  the  bulk  of  the 
electron  density  in  the  phosphorus  lone  pair  is  concentrated  away  from  the  boron 
substituents.  In  this  region,  the  ability  of  the  BH^  group  to  delocalize  the  phosphorus 
lone  pair  is  greatly  diminished.  The  mixing  of  these  orbitals  is  reduced  to  the  point 
that  it  can  be  accurately  examined  using  second  order  perturbation  theory. 

Applying  second-order  perturbation  theory,  the  energy  lowering  that  results 
from  mixing  of  orbitals  is  proportional  to  the  overlap  of  the  orbitals  squared  divided 


81 


by  the  energy  difference  separating  the  orbitals.  It  is  easy  to  show^4  that  the 
numerator  is  twice  as  large  for  HP^H^  as  for  H^PBP^,  which  reflects  the  fact  that 
there  are  two  P-B  interactions  in  PH^BH^.  However,  the  energy  difference  in  the 
denominator  is  not  the  same  in  the  two  molecules,  because  the  in-phase  combination 
of  BHj  p  orbitals  in  HP(BI^)2  is  lower  in  energy  than  the  lone  BP^  p  orbital  in 
HgPBP^.  This  causes  the  denominator  in  the  second-order  perturbation  energy 
expression  to  be  smaller  in  HP^Hj^  than  in  HjPBIij,  so  that  the  energy  lowering 
due  P-B  interactions  are  larger  in  HP(BI^^  than  in  HjPBHj.  Since  in  second-order 
perturbation  theory  interactions  are  additive,  in  the  region  of  very  weak  interactions, 
where  second-order  perturbation  theory  is  valid,  competitive  effects  are  necessarily 
absent. 

Section  3-4-3  Concluding  Remarks  on  the  Plg.^BH^  Molecules 

One  BP^  group  appears  at  first  to  not  have  much  of  an  effect  on  the 
phosphorus  center  since  the  optimized  geometry  is  only  slightly  less  pyramidal  than 
that  of  PHg.  However,  one  BP^  group  actually  has  quite  a  large  effect  on  lowering 
the  barrier  to  planarity  of  the  phosphorus  to  amine-like  values.  A  second  BHj  group 
can  then  form  a  strong  w  bond,  since  the  potential  for  planarizing  phosphorus  has 
already  been  softened  by  the  first.  As  a  result,  the  first  BP^  group  appears  to 
cooperate  with  the  second  by  making  accessible,  without  large  energetic  cost, 
geometries  where  phosphorus  can  form  strong  tt  bonds  to  both  BP^  groups.  Since 
the  presence  of  two  conjugated  BP^  groups  provide  only  a  slightly  softer  potential 


for  planarizing  phosphorus  than  does  one  BF^  group,  addition  of  a  third  BF^  group 
does  not  give  a  cooperative  effect.  Instead,  the  competition  of  the  third  BF^  group 
with  the  first  two  for  electron  density  is  the  dominant  effect;  and  so  the  energy 
required  to  rotate  a  BF^  group  out  of  conjugation  is  smaller  in  P^H^  than  in 
HP(BHV 

These  quantitative  predictions  of  the  sizes  of  the  rotational  barriers,  and  hence 
the  it  bond  strength,  in  Hg.J^BH^  should  be  experimentally  testable  provided  that 
non-sterically  demanding  substituents  could  be  used.  Of  greater  importance  is  the 
qualitative  definition  given  by  this  study,  of  the  conditions  under  which  cooperative, 
rather  than  competitive,  effects  of  electronically  similar,  conjugating  substituents  are 
likely  to  be  observed. 

Section  3-5  Ring  Systems  Containing  P-B  n  Bonds 

In  acyclic  systems,  we  find  that  the  presence  of  one  BHj  group  reduces  the 
barrier  to  planarity  at  phosphorus  but  adding  a  second  BUg  group  nearly  eliminates 
the  barrier.  These  findings  raise  the  question  of  whether  cyclic  (HP-BH\,  systems 
will  have  planar  geometries  at  phosphorus,  since  each  phosphorous  is  flanked  by  two 
borons,  or  whether  they  will  have  non-planar  geometries,  since,  as  the  formula 
implies,  there  is  one  boron  for  each  phosphorous.  Depending  on  whether  n  is  odd 
or  even,  (HP-BH)^  will  have  either  4m  +  2  or  4m  it  electrons  and  hence  will  be, 
respectively,  potentially  aromatic  or  anti-aromatic  in  nature.  This  effect  should  act 
to  make  the  odd  members  of  the  series  tend  to  prefer  geometries  with  planar 
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phosphorus  atoms  and  the  even  members  tend  to  prefer  geometries  with  pyramidal 
phosphorus  atoms. 


Table  3-5.  Boraphosphine  ring  system,  P^H*  and  PjB^He,  relative  using  the 


6-31G*  basis  set. 

Geometrv 

RHF 

MP2 

MP4SDTO 

P2BjH4  Planar  (16) 

30.7 

32.7 

33.3 

PjB^F^  Cis  Hydrogens  (17) 

0 

0 

0 

P^Hj  Trans  Hydrogens  (18) 

6.0 

9.5 

9.2 

P3BsH6  Planar  (19) 

7.4 

0.0 

1.4 

PjRjHg  Chair  (20) 

0 

0 

0 

Indeed,  we  find  for  P^H*  that  a  Qv  geometry,  with  all  hydrogens  cis  and 
pyramidalization  angles  at  phosphorous  of  <p  =  68.2°,  is  preferred  to  a  planar  D2h 
geometry  by  33.3  kcal/mol  at  the  MP4SDTQ  level.  In  contrast,  PsB^Hg  is  computed 
to  have  a  smaller  pyramidalization  angle,  <p  =  59.0°  at  the  RHF  level,  and  the 
barrier  to  planarization  at  all  three  phosphorus  atoms  (Qv  D3h)  is  only  1.4 
kcal/mol  at  the  MP4SDTQ  level. 

The  aromatic  system  has  such  a  low  barrier  to  planarization  that  even  better 
calculations  might  find  a  planar  equilibrium  geometry.  As  noted  in  Section  3-1-3,  a 
planar  equilibrium  geometry  has  been  found  for  a  highly  substituted  derivative  of 
(IIPBH)^  by  x-ray  crystallography.19  Since  our  calculations  were  performed  on  the 
parent  molecule  and  refer  to  the  gas-phase,  the  discrepancy  between  our 
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computational  result  and  experiment  could,  of  course,  be  due  to  substituent  effects 
or  to  packing  forces  in  the  crystal. 

Although  calculations  predict  the  Qv  structure  (17)  of  PjB^  to  be  the 
ground  state,  X-ray  crystallography  of  a  substituted  derivative  resembles  Qh  structure 
(18),  which  we  calculate  to  be  9.2  kcal/mol  higher  in  energy  at  the  MP4SDTQ  level. 
Examination  of  structure  17  reveals  that  all  the  substituents  would  be  confined  to  the 
same  side  of  the  PB  ring.  Since  the  substituted  derivative  contains  sterical)'' 
demanding  groups,20  this  may  account  for  the  fact  that  the  less  sterically  crowded  18 
is  found  to  be  preferred  for  this  derivative.  However,  we  expect  a  derivative  with 
less  sterically  demanding  substituents  to  adopt  a  conformation  similar  to  structure  17. 
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Chapter  4  Study  of  Oxyallyl  and  Dimethyl  Oxyallyl 


Section  4-1  Introduction 

Oxyallyl  and  dimethyloxyallyl  are  molecules  in  which  two  electrons  occupy  two 
orbitals  that  are  close  in  energy.  Thus,  these  are  two  examples  of  diradicals.1  Since 
no  classical  Kekul6  structures  can  be  written  for  these  diradicals,  they  are  also  called 
non-Kekul6  molecules.2,3,4 

Hund’s  rule5  predicts  that  the  lowest  energy  configuration  for  two  electrons 
occupying  two  degenerate  orbitals  will  be  that  which  produces  the  greatest  spin 
multiplicity.  This  rule  was  first  applied  to  atoms,  but  it  is  also  valid  for  many 
molecules.1  However,  as  discussed  by  Borden1  and  Du, 3  diradicals  in  which  the 
unpaired  electrons  are  localized  to  different  regions  in  space  can  have  a  singlet 
ground-state  wave  function,  even  though  Hund’s  rule  predicts  that  the  triplet  should 
be  the  ground  state.  The  reason  why  such  diradicals  might  be  expected  to  violate 
Hund’s  rule  is  that  the  Coulomb  repulsion  energy  of  the  two  antiparallel  spin 
electrons  is  significantly  reduced  if  they  are  each  confined  to  a  different  region  of 
space.  Electron  correlation  between  these  two  electrons  and  the  rest  of  the  electrons 
then  becomes  the  determining  factor  in  whether  the  singlet  or  triplet  wave  function 
is  lower  in  energy;  and  the  singlet  is  usually  favored. 

As  a  result,  in  calculations  of  singlet-triplet  energy  separations  in  diradicals, 
electron  correlation  is  very  important.  Additionally,  since  singlet  diradicals  cannot 
be  adequately  described  by  a  single  Slater  determinant,  Mtfllcr-Plesset  perturbation 
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theory  will  not  provide  accurate  results.  Consequently,  MCSCF  or  Cl  calculations 
are  required. 

Section  4-1-2  Prior  Work  on  Oxyallyl 

Oxyallyl  (2)  is  formed  as  an  intermediate  in  the  reaction  that  takes  allene 
oxides  to  cyclopropanone^  and  in  the  isomerization  and  cycloaddition  reactions  of 
cyclopropanone7.  Since  oxyallyl  is  isoelectronic  with  trimethylenemethane  (TMM) 
(1),  both  might  be  expected  to  have  a  ground  state  of  the  same  spin  multiplicity. 
Experiments  have  found  that  TMM  and  derivatives  have  a  triplet  ground  state,8,9 
which  agrees  with  theoretical  predictions.1  However,  we  are  unaware  of  any 
experimental  evidence  of  a  triplet  ground-state  for  oxyallyl  or  derivatives. 


Figure  4-1  Structure  of  TMM,  Oxyallyl  and  Dimethyloxyallyl. 
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Several  ab  initio  studies  have  predicted  that  oxyallyl  has  a  triplet  ground 
state.10,11  MCSCF  calculations  found  that  the  singlet-triplet  gap  to  be  12.3  kcal/mol 
with  a  3-2 1G  basis  set  and  5.6  kcal/mol10  when  polarization  functions  were  added. 
The  large  change  in  the  energy  gap  upon  adding  polarization  functions  suggested  to 
the  authors  that,  had  polarization  functions  been  present  during  geometry 
optimization,  the  singlet  might  have  been  stabilized  even  more,  relative  to  the  triplet; 
but  the  authors  did  not  explore  this  possibility  computationally.  Apparently,  the 1  Aj 
state  has  a  strong,  polar  carbonyl  bond  that  requires  polarization  function  for  an 
adequate  description.12 

Recent  theoretical  work  on  2,4-dimethylenecyclobutane- 1 ,3  diyl  and  the  mono- 
and  dioxo  derivatives13  suggests  that  as  in  1,  substitution  of  O  for  CTLj  selectively 
stabilizes  the  singlet.  In  fact,  the  dioxo  derivative  was  predicted  to  have  a  singlet 
ground  state.  This  study  found  that  in  all  three  diradicals,  correlation  of  the  a  with 
the  n  electrons  selectively  stabilizes  the  singlet  state. 

The  addition  of  a  -  n  electron  correlation  and  inclusion  of  polarization 
functions  during  geometry  optimization  should  reduce  the  theoretically  predicted 
singlet-triplet  gap  in  oxyallyl,  possibly  leading  to  the  finding  that  oxyallyl  has  a  ground 
state  singlet.  We  therefore  attempted  to  recalculate  the  ground  state  of  oxyallyl, 
using  a  better  basis  set  and  including  more  extensive  electron  correlation. 
Dimethyloxyallyl  was  also  chosen  for  study  since  the  JAj  state  of  oxyallyl  appears  to 
have  a  contribution  from  a  zwitterionic  resonance  structure,  2b.  If  this  is  true,  then 


addition  of  methyl  groups  to  the  terminal  carbons  should  selectively  stabilize  the 
singlet,  relative  to  the  triplet. 

Section  4-2  Computational  Methodology 

Calculations  were  carried  out  using  the  Dunning  double  zeta  basis  set14  with 
polarization  functions  on  all  non-hydrogen  atoms.16  Complete  Active  Space  (CAS) 
MCSCF  calculations,  using  all  the  spin  adapted  configurations  (12  for 1  \  and  9  for 
8B^)  that  arise  from  4  ir  electron  occupying  the  lowest  4  jt  orbitals,  were  performed 
to  optimize  the  geometry  of  oxyallyl.  The  MCSCF  calculations  were  performed  at 
the  Institute  for  Molecular  Science  in  Okazaki,  Japan  using  the  MOLPRO  package 
of  ab  initio  programs.18  Contracted  CI-SD  calculations17  were  also  performed  with 
MOLPRO.  Strictly  variational  Cl  calculations  were  carried  out  using  the  same  DZP 
basis  set  for  carbon  and  oxygen  and  an  improved  triple  zeta  basis  set  for 
hydrogen.18  These  Cl  calculations  were  performed  at  the  University  of  Washington 
using  the  MELD  package  of  ab  initio  programs.19 

The  C-H  bond  lengths  and  bond  angles  were  fixed  at  the  values  optimized  in 
the  earlier  MCSCF  study.10  The  C-O  and  C-C  bond  lengths  were  optimized  by  a 
quadratic  2-dimensional  fit  of  the  CAS-MCSCF  energies,  using  6  points  for  each 
state.20 

For  dimethyl  oxyallyl,  starting  from  the  MCSCF  optimized  oxyallyl  geometry, 
the  hydrogen  cis  to  the  oxygen  was  replaced  with  methyl  groups  at  a  fixed  HgC-C 
bond  length  of  1.499 A.  The  other  geometric  parameters  for  the  methyl  groups  were 
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MCSCF  Optimized  geometry  of  Oxyallyl. 
'At  State  and  3B2  State 

O 


Figure  4-2  Geometric  data  for  Oxyallyl 

taken  from  the  results  of  an  UHF/6-31G*  calculations  on  the  ethyl  radical*1  The 
methyl  group  conformation  was  such  that  the  hydrogen  in  the  molecular  plane 
pointed  toward  the  oxygen.  A  ROHF  triplet  calculation  showed  this  conformation 
was  lower  in  energy  by  0.0004  hartree  (0.25  kcal/mol)  than  one  with  in  which  the 


hydrogen  pointed  away  from  the  oxygen. 
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Section  4-2-1  Cl  Configuration  Selection 

The  Cl  calculations  on  both  oxyallyl  and  dimethyloxyallyl  started  with 
molecular  orbitals  derived  from  a  RHF  wave  function  for  the  triplet  and  a  TCSCF 
wave  function  for  the  singlet.  The  virtual  orbitals  were  then  transformed  into  K 
orbitals,22  in  order  to  generate  virtual  orbitals  that  would  provide  the  maximum 
amount  of  correlation  with  the  fewest  numbers  of  orbitals.  All  Is  orbitals  on  carbon 


Table  4-1.  Number  of  spin  adapted  configurations  considered  in  Cl  calculations 
on  oxyallyl  and  dimethyloxyallyl. 


Number  of  Spin  Adapted  Configurations 


Cl  Level 

Oxyallyl 

Dimethyloxyallyl 

% 

% 

‘*1 

TCSCF/RHF 

2 

1 

2 

1 

tt-CAS  “ 

12 

9 

900 

1192 

Full  7T 

6708 

9024 

tt-CAS  +  tt-SD 

1234 

1580 

350472 

592325 

tt-CAS  +  a-Sb 

8679 

12645 

345176 

600942 

tt-CAS +  a-S, tt-SD 

9901 

14271 

“CAS  for  Oxyallyl  is  all  configurations  of  the  4  n  electrons  in  the  4  lowest  n 
orbitals.  For  Dimethyloxyallyl,  all  configurations  of  the  8  *  electrons  in  the  8 
lowest  ir  orbitals. 

bFor  Dimethyloxyallyl,  a  excitations  into  the  a  virtual  orbitals  have  been  limited 
to  only  the  lowest  13  unoccupied  K  orbitals. 
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and  oxygen  were  cored.  Several  different  types  of  Cl  calculations  were  then 
performed,  starting  with  a  complete  active  tt  space  (rr-CAS)  Cl  calculation  for  4  n 
electrons  in  4  tt  orbitals  for  oxyallyl  and  8  jr  electrons  in  8  n  orbitals  for 
dimethyloxyallyl.23  The  number  of  configurations  for  each  Cl  calculation  appears 
in  Table  4-1  and  the  absolute  energies  of  the  1Ai  states  are  reported  in  Table  4-2. 

Table  4-2.  Energies  of  1A1  states  of  oxyallyl  and  dimethyl  oxyallyl  computed  at 
various  levels  of  theory  expressed  in  hartrees. 


Cl  Level 

Oxyallyl 

Dimethyloxyallyl 

TCSCF/RHF 

-190.7464 

-268.8156 

7T-CASa 

-190.7684 

-268.8775 

Full  jr 

-190.7886 

tt-CAS  +  n- SD 

-190.7886 

-268.9067 

ir-CAS  +  a-Sb 

-190.8677 

-268.9721 

jt-CAS +  a-S,7r-SD 

-190.8798 

aCAS  for  Oxyallyl  is  all  configurations  of  the  4  jr  electrons  in  the  4  lowest  n 
orbitals.  For  Dimethyloxyallyl,  all  configurations  of  the  8  rr  electrons  in  the  8 
lowest  k  orbitals. 

bFor  Dimethyloxyallyl,  a  excitations  into  the  a  virtual  orbitals  have  been  limited 
to  only  the  lowest  13  unoccupied  K  orbitals. 

Additional  jr  correlation  was  added  to  the  jt-CAS  Cl  wave  function  by 
inclusion  of  all  single  and  double  excitations  in  the  jr  space  from  all  the  reference 
jt-CAS  configurations  (jt-CAS  +  jt-SD).  To  verify  the  jt-CAS  +  jt-SD  calculation 
recovered  nearly  all  of  the  possible  correlation  energy  in  the  ir  space,  FULL-jt  Cl 
was  performed  on  oxyallyl.  Next,  correlation  between  the  a  and  jt  electrons  was 
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added  by  including  single  a  excitations  from  all  7T-CAS  configurations  (ir-CAS  + 
a-S).  Because  of  the  size  of  dimethyloxyallyl,  the  a-S  Cl  calculation  had  to  be 
limited  to  only  excitations  into  the  lowest  13  a  virtual  K  orbitals.  To  get  a  measure 
of  the  interaction  between  n  and  a-ir  correlation,  an  additional  Cl  was  performed  on 
oxyallyl  that  allowed  either  single  excitation  within  the  a  space  or  single  and  double 
excitations  within  the  it  space  from  all  rr-CAS  configurations  (jt-CAS  +  a-S,  u-SD). 

For  oxyallyl,  it  was  also  possible  to  perform  SD-CI,  using  the  MCSCF 
molecular  orbitals  and  allowing  all  possible  a  and  n  single  and  double  excitations 
from  the  tt-CAS  reference  space.  The  energies  of  these  Cl  wave  functions  were 
calculated  with  the  internally  contracted  Cl  method  of  Werner  and  Knowles,16  which 
assumes  certain  relationships  between  the  coefficients  of  sets  of  configurations  and 
thus  reduces  the  size  of  the  Cl  calculation  that  is  actually  performed.  The  MCSCF 
and  contracted  Cl  energies  are  reported  in  Table  4-3.  However,  it  should  be  noted 
that,  because  the  relative  weights  of  the  configurations  in  some  groups  are  assumed, 
the  contracted  Cl  calculations  are  not  strictly  variational. 

Section  4-3  Results 
Section  4-3-1  Oxyallyl 

The  MCSCF  geometries  provided  by  reoptimizing  the  C-O  and  C-C  bond 
lengths  with  using  the  DZP  basis  set  resulted  in  a  significantly  shorter  C-O  bond 
length  than  when  the  3-21G  basis  set  was  employed.  The  C-O  bond  in  the  8Bj  state 
shortened  by  nearly  0.059  A  with  an  accompanying  increase  in  each  of  the  C-C  bonds 
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by  half  this  amount,  0.029  A.  For  the  state,  only  the  C-O  bond  length  changed; 
it  shortened  by  0.020  A.  The  C-O  bond  shortening  is  not  surprising  since  adding 
polarization  functions  should  describe  the  carbonyl  bona  better. 

Using  the  original  3-21G  optimized  geometries,  the  CAS-n  MCSCF  energy 
difference  of  2.8  kcal/mol  between  the  singlet  and  the  triplet  with  the  Dunning  DZP 
basis  set  is  only  half  of  the  5.6  kcal/mol  computed  with  the  3-21 G*  basis  set.10  Since 
geometry  reoptimization  affects  the  3Rj  geometry  more  than  the  1A1,  the  MCSCF 
value  of  A  Es_T  increases  to  4.4  kcal/mol  upon  geometry  reoptimization. 


Table  4-3.  Results  of  MCSCF  and  contracted  Cl  calculations  on  oxyallyl. 


Number  of  spin  adapted 
configurations 


Calculation* 

aA, 

% 

Energy  of  Aj 
in  hartrees 

A  fcg_*p  in 
kcal/mol 

MCSCF 

12 

9 

-190.7580 

4.4 

CI-SD 

812362 

1480927 

-190.2072 

2.4 

CI-SDQ^ 

-190.2703 

1.6 

aAll  calculations  were  carried  out  using  a  CAS-tt  MCSCF  wave  function  as  the 
reference  space.  All  single  and  double  excitations  were  added  using  the 
contracted  Cl  method  of  Werner  and  Knowles. 

bQuadruple  excitations  were  added  using  the  Davidson  method. 


The  MCSCF  value  of  4.4  kcal/mol  for  AE^_Tis  slightly  larger  than  the  2.3 
kcal/mol  obtained  from  the  CAS-rr  Cl  calculations.24  The  CAS-w  Cl  calculations, 
which  used  the  same  set  of  configurations  as  the  MCSCF,  do  not  fully  regain  as  much 
correlation  energy  since  the  virtual  orbitals  have  not  been  optimized  variationally  to 
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provide  the  largest  amount  of  ir  electron  correlation.  Since  the  triplet  wave  function 
is  more  delocalized  than  the  singlet,25  as  expected,15  improvements  in  correlation 
of  the  n  electrons  selectively  stabilize  the  triplet,  relative  to  the  singlet.  In  fact,  since 
the  RHF/TCSCF  calculation  provides  correlation  for  only  the  two  non-bonding  n 
electrons,  it  actually  predicts  a  singlet  ground  state.  As  shown  in  Table  4-3  and 
Table  4-4  respectively,  addition  of  ir  correlation  of  the  non-bonding  and  n  bonding 
electrons  at  either  the  CAS-tt  MCSCF  or  tt-CAS  +  tt-SD  Cl  level  selectively 
stabilizes  the  triplet  by  9.5  kcal/mol  as  compared  to  the  RHF/TCSCF  energy. 


Table  4-4.  Predicted  AE^.Tfor  oxyallyl  and  dimethyloxyallyl  at  various  levels  of 
Cl  theory. 


AE^.t  (kcal/mol) 


Cl  Level 

Oxyallyl 

Dimethyloxyallyl 

TCSCF/RHF 

-5.1 

-9.8 

tt-CAS* 

2.3 

00 

ri 

i 

Full  ir 

4.4 

tt-CAS  +  ir- SD 

4.4 

-1.3 

tt-CAS  +  cr-Sb 

-0.3 

-5.3 

tt-CAS  +  a-S,  tt-SD 

0.8 

“CAS  for  Oxyallyl  is  all  configurations  of  the  4  ir  electrons  in  the  4  lowest  ir 
orbitals.  For  Dimethyloxyallyl,  all  configurations  of  the  8  ir  electrons  in  the  8 
lowest  ir  orbitals. 

bFor  Dimethyloxyallyl,  a  excitations  into  the  a  virtual  orbitals  have  been  limited 
to  only  the  lowest  13  unoccupied  K  orbitals. 
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Also  as  expected,13  the  more  ionic  n  wave  function  for  the  singlet  state  is 
selectively  stabilized  through  the  inclusion  of  a-n  correlation.  Inclusion  of  a-S  in 
either  the  ir-C AS  Cl  or  the  tt-CAS  +  n-SD  Cl  calculations  predicts  AEg.jto  be  near 
zero.  The  tt-CAS  +  jt-SD,  a-S  value  of  0.8  kcal/mol  is  in  good  agreement  with  the 
value  in  Table  4-3  of  1.6  kcal/mol  for  the  tt-CAS  +  SD-CI  calculation  obtained  with 
contracted  Cl  methodology.  The  exact  value  of  AEg_T  in  oxyallyl  is  less  important 
than  the  fact  that  it  is  almost  an  order  of  magnitude  smaller  than  AEg_x  in  the 
isoelectric  TMM  (l).1 

An  inspection  of  the  n  population  analyses26  for  2  in  Table  4-5  supports  the 
assertion  that  the  l\  has  more  charge  separation  than’Bj.  Contributions  from  the 
zwitterionic  structure  2b,  could  be  responsible  for  the  reduction  in  AEs_T  in  2, 
compared  to  1.  However,  although  the  greater  electron  density  on  oxygen  in  lAlf 
compared  to  the  3Bj  supports  this  explanation,  the  shorter  C-O  and  longer  C-C  bond 
lengths  in  2Aj  do  not.  If  structure  2b  was  very  important  in  1Al,  one  would  expect 
the  C-0  bond  length  to  be  longer  and  the  C-C  bond  lengths  to  be  shorter  than  in  the 

3b2. 

Earlier  work  suggested  the  major  difference  between  1  and  2, 10  is  the  greater 
strength  of  a  C-O  v  bond,  relative  to  a  C-C  *  bond.  On  going  from  1  to  2,  the 1  Ax 
state  is  stabilized  relative  to  the  other  low  lying  singlet  states,27  2Bj  and  1Bt,  because 
the  *Aj  wave  function  is  the  only  one  that  allows  a  C-O  *  bond  to  form.  Since  the 
3Bj  wave  function  could  be  thought  of  as  a  mixture  of  all  three  possible  covalent 
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resonance  structures,  only  one  of  which  has  a  full  C-O  tt  bond,  it  is  stabilized  less 
than  1A1  state  by  the  substitution  of  O  for  CH^. 

Table  4-5.  tt  atomic  orbital  total  and  (unpaired)  electron  population^’1*  in  the 
and  3Bj  states  of  oxyallyl  at  various  levels  of  theory. 


Calculation 

State 

o 

Q(Q) 

C2 

TCSCF 

1.42 

0.93 

0.68 

RHF 

3B> 

1.33  (0.22) 

0.94  (0.87) 

0.77  (0.03) 

CAS-tt  +  tt-SD 

1.39 

0.91 

0.77 

CAS-tt  +  tt-SD 

% 

1.21  (0.48) 

0.96  (0.85) 

0.86  (-0.18) 

CAS-tt  +  a-S,  tt-SD 

l\ 

1.38 

0.89 

0.82 

CAS-tt  +  a-S,  tt-SD 

3B2 

1.20  (0.45) 

0.96  (0.83) 

0.87  (-0.13) 

aSee  reference  27. 


bBecause  there  is  some  electron  density  in  d-*r  orbitals,  the  p-jr  total  populations 
do  not  sum  to  4.0;  nor  do  the  unpaired  electron  populations  in  the’B^  state  sum 
exactly  to  2.0. 


Another  difference  between  1  and  2  is  that  symmetry  mandates  that  all  three 
resonance  structures  contribute  equally  to  the  triplet  state  in  1,  but  not  in  2.  Indeed, 
as  shown  in  Table  4-5,  the  ratio  of  unpaired  spin  density  on  carbon  to  that  on  oxygen 
ranges  from  1.77  with  tt-CAS  +  tt-SD  to  1.84  with  tt-CAS  +  a-S,  tt-SD  Cl.  This 
ratio  differs  significantly  from  the  ratio  of  1.00  demanded  by  symmetry  in  1.  From 
the  ratio  in  2,  it  can  be  calculated  that  2a  contributes  about  2.6  times  more  than 
either  of  the  two  resonance  structures  that  containing  a  C-C  tt  bond. 

One  way  to  understand  the  near  degeneracy  between  the’B,  and  1A1  states 
of  oxyallyl  begins  with  the  fact  that  structure  2a  represents  the  bulk  of  the  bonding 
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in  both  states.  The  triplet  does  have  two  other  covalent  structures  that  are  minor 
contributors,  and  the  more  delocalized  wave  function  for  the  triplet  account  for  the 
differences  in  optimized  bond  length  between  the  two  states.  However,  lAj  has  a 
contribution  from  2b,  which  is  not  a  contributor  to  the  triplet;  and  this  accounts  for 
the  differences  in  the  tt  electron  densities  between  the  two  states. 

The  triplet  is  selectively  stabilized  by  a  more  delocalized  n  wave  function. 
The  singlet  is  selectively  stabilized  oxygen’s  ability  to  accept  excess  electron  density. 
These  two  effects  are  apparently  of  about  the  same  magnitude,  which  results  in  near 
degeneracy  of  the  energies  of  the3^  and  1\  in  oxyallyl. 

Section  4-3-2  Dimethyloxyallyl 

Because  singlet  oxyallyl  is  selectively  stabilized  through  the  contribution  of  2b, 
addition  of  methyls  might  result  in  a  singlet  ground  state.  This  is  confirmed  by  Cl 
calculations  on  dimethyloxyallyl  (3),  the  results  of  which  are  reported  in  Table  4-4. 
At  all  levels  of  Cl  theory,  AE%_T  for  3  is  smaller  than  for  2.  The  negative  sign  for 
every  value  of  AE%_Tin  3  indicates  that,  at  all  levels  of  theory,  a  singlet  ground  state 
is  predicted. 

Comparison  of  the  A  E^_T  values  at  the  same  levels  of  theory  for  2  and  3  show 
a  5  to  6  kcal/mol  stabilization  of  1A1,  relative  to 3 B^,  in  3.  This  finding  suggests  that 
if  a  jr-CAS  +  ct-S,  jt-SD  or  SDQ  Cl  calculation  were  performed  on  3,  the  AEs„Tfor 
3  would  be  around  -4  to  -5  kcal/mol.  Indeed,  comparisons  of  Cl  calculations  with 
partial  and  complete  a  single  excitations  from  a  smaller  set  of  reference 
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configurations  indicate  that  had  all  a  virtual  orbitals  been  used  for  correlation  in  the 
tt-CAS  +  ct-S  Cl  calculation,  the  singlet  would  have  been  selectively  stabilized  by 
another  1  to  2  kcal/mol.  Therefore,  the  5  kcal/mol  appears  to  be  a  conservative 
estimate  of  the  magnitude  by  which  1Al  lies  below  3B^  in  3. 

Section  4-4  Conclusion 

Although  our  calculations  on  oxyallyl  predict  a  small  preference  for  the  triplet, 
the  calculations  on  dimethyloxyallyl  unequivocally  predict  a  singlet  ground  state.  This 
result  is  consistent  with  predictions  of  Lahti  and  coworkers?8  on  cyclic  derivatives 
of  oxyallyl.  Our  results  indicate  that  the  constraints  on  the  C-C-C  bond  angle  at  the 
carbonyl  group,  imposed  by  incorporation  of  oxyallyl  in  a  four-  or  five  member  ring, 
are  not  required  to  obtain  a  singlet  ground  state  for  dialkyl  derivatives  of  oxyallyl.29 
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APPENDIX  1.  Geometries  and  Energies  of  Carbon  Centered  Radicals. 
Bond  Lengths  in  A,  Angles  in  Degrees,  Energies  in  Hartrees. 


Substituent,  X 

C-Hl 

C-H2 

C-X 

X-H3 

X-H4 

H3-X-H4 

EOJHF/6-31G* 

1X-C-H1 

X-C-H2 

H-C-H 

C-X-H3 

C-X-H4  H3-X-H5 

Li 

1.092 

1.092 

1.943 

-46.39468 

126.6 

126.6 

106.8 

BeH 

1.082 

1.082 

1.665 

1.336 

. 

-54.19425 

124.6 

124.6 

110.8 

180.0 

BH2  rot 

1.079 

1.079 

1.555 

1.193 

1.193 

117.2 

-64.80976 

122.5 

122.5 

115.0 

121.4 

121.4 

bh2 

1.078 

1.078 

1.536 

1.190 

1.190 

119.0 

-64.82416 

123.0 

123.0 

114.0 

120.5 

120.5 

CH, 

1.075 

1.075 

1.498 

1.086 

1.086 

108.0 

-78.59715 

120.4 

120.4 

117.3 

111.4 

111.4 

107.1 

118 


Substituent,  X 

C-Hl 

C-H2 

C-X 

X-H3 

X-H4 

H3-X-H4 

X-C-H2 

H-C-H 

C-X-H3 

C-X-H4  H3-X-H5 

NHg  rot 

1.079 

1.074 

1.420 

1.002 

1.003 

107.1 

-94.55750 

120.3 

120.3 

119.2 

111.4 

111.4 

NH, 

1.076 

1.076 

1.402 

0.999 

0.999 

109.5 

-94.58673 

115.7 

115.7 

116.9 

113.5 

113.5 

NHg-t 

1.070 

1.070 

1.470 

1.013 

1.013 

107.2 

115.7 

115.7 

124.3 

107.2 

107.2 

105.8 

OH  rot 

1.076 

1.076 

1.366 

0.948 

-114.40433 

116.1 

116.1 

119.2 

110.5 

OH 

1.078 

1.073 

1.359 

0.946 

60.2 

-114.40876 

112.8 

117.6 

118.6 

110.2 

NC 

1.071 

1.071 

1.347 

1.167 

-131.27045 

118.7 

118.7 

112.7 

180.0 

CN 

1.072 

1.072 

1.391 

1.159 

-131.30689 

119.9 

119.9 

120.2 

180.0 

F 

1.073 

1.073 

1.331 

-138.40211 

113.9 

113.9 

121.3 

CHO 

1.072 

1.074 

1.414 

1.237 

1.086 

119.3 

-152.30039 

119.5 

121.2 

119.9 

122.2 

120.0 

-152.30039 


119 


Substituent,  X 

C-Hl 

C-H2 

C-X 

X-H3 

X-H4 

H3-X-H4 

1X-C-H1 

X-C-H2 

H-C-H 

C-X-H3 

C-X-H4  H3-X-H5 

NO, 

1.068 

1.068 

1.418 

1.198 

1.198 

125.8 

-243.02982 

116.5 

116.5 

126.9 

117.1 

117.1 

SiH, 

1.078 

1.078 

1.860 

1.477 

1.477 

109.5 

-329.64374 

122.6 

122.6 

114.7 

110.2 

110.2 

107.9 

PH,  rot 

1.075 

1.075 

1.821 

1.405 

1.405 

95.3 

-380.85599 

121.5 

121.5 

117.0 

98.8 

98.8 

PH, 

1.076 

1.076 

1.807 

1.403 

1.403 

96.2 

-380.85710 

119.9 

119.9 

116.7 

99.4 

99.4 

SH  rot 

1.086 

1.086 

1.715 

1.373 

-436.95463 

109.6 

109.6 

1103 

117.2 

SH 

1.073 

1.072 

1.749 

1.326 

-437.07234 

116.0 

120.1 

119.4 

Cl 

1.070 

1.070 

1.718 

-498.46108 

116.5 

116.5 

122.3 

Appendix  2.  Geometries  and  Energies  of  Silicon  Centered  Radicals. 


Bond  Lengths  in  A,  Angles  in  Degrees,  Energies  in  Hartrees. 


Substituent,  X 

SiHl 

SiH2 

SiX 

X-H3 

X-H4 

H3-X-H4 

E(UHF/6-31G*l 

XSIHl 

XSiH2 

HSiH 

SiXH3 

SiXH4  H3-X-H5 

Li 

1.495 

1.495 

2.426 

-297.47969 

128.1 

128.1 

102.9 

BeH 

1.497 

1.497 

2.166 

1334 

-305.23926 

120.2 

120.2 

108.8 

180.0 

BH,  rot 

1.469 

1.469 

2.045 

1.188 

1.188 

117.5 

-315.83391 

121.8 

121.8 

116.4 

1213 

121.3 

BH,  rot  A  pyr 

1.482 

1.482 

2.046 

1.189 

1.191 

117.3 

-315.84524 

1113 

111.5 

108.7 

1223 

120.2 

BH, 

1.471 

1.471 

1.954 

1.188 

1.188 

119.6 

-315.85687 

124.0 

124.0 

112.0 

120.2 

120.2 

CH, 

1.480 

1.480 

1.894 

1.085 

1.085 

107.8 

-329.65181 

111.6 

111.6 

109.2 

110.9 

110.9 

107.9 

Substituent,  X 

SiHl 

SiH2 

SiX 

X-H3 

X-H4 

H3-X-H4 

i  XSIH1 

X§iH2 

HSiH 

SiXH3 

SiXH4 

H3-X-H5 

NH2  rot 

1.482 

1.721 

0.995 

0.998 

112.5 

-345.66023 

113.3 

113.3 

107.7 

123.2 

123.2 

NHj 

1.480 

1.729 

0.998 

0.998 

110.9 

-345.66304 

108.5 

110.2 

120.2 

120.2 

NHj  + 

1.461 

1.461 

1.928 

1.012 

1.012 

106.7 

-346.01412 

101.9 

116.5 

112.2 

112.2 

106.6 

OH  rot 

1.476 

1.476 

1.652 

0.945 

-365.50603 

110.4 

109.3 

118.7 

OH 

1.483 

1.473 

1.652 

0.946 

59.8 

-365.50694 

111.9 

106.9 

109.9 

118.2 

NC 

1.470 

1.740 

1.163 

-382.35293 

108.3 

112.1 

180.0 

CN 

1.468 

1.468 

1.853 

1.142 

-382.35497 

108.6 

113.4 

180.0 

F 

1.475 

1.475 

1.599 

-389.52283 

108.8 

110.5 

CHO 


1.462  1 

120.5  1 


-403.31552 


.464  1.863 

19.2  120.3 


1.215 

122.4 


1.093  119.5 

118.1 


122 


Substituent,  X 

SiHl 

SiH2 

SiX 

X-H3 

X-H4 

H3-X-H4 

FiTJHF/6-31G*'l 

i  XSIH1 

XSiH2 

HSiH 

SiXH3 

SiXH4  H3-X-H5 

NQj 

1.464 

1.464 

1.873 

1.197 

1.197 

125.1 

-494.08775 

106.0 

106.0 

114.9 

117.5 

117.5 

SiH, 

1.479 

1.479 

2.345 

1.477 

1.477 

108.8 

-580.69070 

113.6 

113.6 

109.4 

109.4 

109.4 

108.5 

PH,  rot 

1.478 

1.475 

2.271 

1.404 

1.403 

94.5 

-631.91699 

112.4 

112.4 

109.5 

97.3 

97.3 

PH, 

1.478 

1.478 

2.264 

1.403 

1.403 

95.8 

-631.91914 

109.2 

109.2 

109.2 

98.0 

98.0 

SH  rot 

1.474 

1.474 

2.156 

1.328 

-688.15082 

110.8 

110.8 

109.6 

97.9 

SH 

1.477 

1.473 

2.149 

1.328 

-688.15287 

112.0 

105.6 

110.6 

97.9 

Cl 

1.472 

1.472 

2.070 

-749.56312 

108.7 

108.7 

111.2 

APPENDIX  3.  Geometries  and  Energies  of  Selected  Closed-Shell  Carbon  Molecules* 
Bond  Lengths  in  A,  Angles  in  Degrees,  Energies  in  Hartrees. 


Substituent,  X 
EflRHF/6-3lG' 

C-Hl 

’1X-C-H1 

C-H2 

X-C-H2 

c-x 

H-C-H 

X-H3 

C-X-H3 

X-H4  H3-X-H4 

C-X-H4  H3-X-H5 

CN 

1.082 

1.082 

1.468 

1.135 

-141.92753 

109.8 

109.8 

1093 

180.0 

NC 

1.081 

1.081 

1.421 

1.153 

-131.89436 

109.6 

109.6 

1093 

180.0 

NQ, 

1.080 

1.076 

1.478 

1.191 

1.191 

125.1 

-243.66199 

106.5 

108.0 

112.9 

108.0 

108.0 

iJh, 

1.078 

1.078 

1.507 

1.012 

1.012 

107.3 

-95.90188 

108.1 

108.1 

110.8 

1166 

116.6 

107.3 

‘Numbers  not  listed  appear  in  the  Carnegie-Mellon  Quantum  Chemistry  Archive  or 
in  Luke,  B.  T.;  Pople,  J.  A;  Krogh-Jespersen,  M.-B.  Apeloig,  Y.;  Chandrasekhar,  J.; 
Schleyer,  P.  v.  R.  J.  Am.  Chem.  Soc.  1986, 108 ,  260. 


APPENDIX  4.  Geometries  and  Energies  of  Selected  Closed-Shell  Silicon  Molecules* 
Bond  Lengths  in  A,  Angles  in  Degrees,  Energies  in  Hartrees. 


Substituent,  X 

C-Hl 

C-H2 

C-X 

X-H3 

X-H4 

H3-X-H4 

FfRHF/6-3lG 

*1X-C-H1 

X-C-H2 

H-C-H 

C-X-H3 

C-X-H4  H3-X-H5 

CN 

1.467 

1.467 

1.866 

1.138 

-382.97440 

107.5 

107.5 

1113 

180.0 

NC 

1.467 

1.467 

1.742 

1.162 

-382.97508 

108.1 

108.1 

110.8 

180.0 

NO, 

1.463 

1.461 

1.870 

1.197 

1.197 

127.0 

-494.70911 

103.7 

106.1 

112.6 

1173 

117.5 

CHO 

1.469 

1.475 

1.927 

1.193 

1.100 

118.4 

-4C3.93843 

109.5 

108.2 

108.9 

122.7 

118.9 

‘Numbers  not  listed  appear  in  the  Carnegie-Mellon  Quantum  Chemistry  Archive  or 
in  Luke,  B.  T.;  Pople,  J.  A*;  Krogh-Jespersen,  M.-B.  Apeloig,  Y.;  Chandrasekhar,  J.; 
Schleyer,  P.  v.  R.  /.  Am.  Chem.  Soc.  1986, 108 ,  260. 


Appendix  5.  Geometries  and  Energies  of  Saturated  Allylic  Systems 


H 

I  xH? 

H2^X^z>2 

A,  A, 


Bond  Lengths  are  in  angstroms  and  angles  in  degrees. 

X  Y  Z 

XY  YZ 

H1X  H2X 

HY 

H1Z 

H2Z 

XYZ 

H1XYH2XY 

HYX 

H1ZY 

H2ZY 

ECRHF/6-31G*) 

HXH 

HYZ 

H1ZH2 

H2ZH2 

_  _  _ 

-  - 

-  - 

— 

— 

— 

c  c  c 

1.319  1.503 

1.077  1.076 

1.079 

1.084 

1.087 

125.3 

121.9  121.7 

118.9 

111.4 

110.9 

-117.071471 

116.5 

115.8 

108.2 

107.0 

C  C  Si 

1.325  1.873 

1.077  1.078 

1.081 

1.477 

1.478 

123.7 

122.3  122.3 

117.7 

109.2 

110.9 

-368.112508 

115.5 

118.6 

109.1 

107.6 

C  Si  C 

1.693  1.879 

1.077  1.076 

1.471 

1.085 

1.086 

125.4 

123.2  121.9 

120.6 

110.6 

111.0 

-368.085435 

114.9 

114.0 

108.2 

107.6 

Si  C  C 

1.700  1.508 

1.468  1.467 

1.080 

1.085 

1.088 

128.2 

123.0  121.9 

116.9 

112.3 

111.0 

-368.071536 

115.1 

114.8 

107.7 

107.0 

Si  Si  C 

2.127  1.887 

1.469  1.468 

1.471 

1.084 

1.085 

125.6 

122.6  121.9 

119.6 

110.7 

110.8 

-619.122356 

115.5 

1T4.9 

108.2 

108.0 

Si  C  Si 

1.696  1.855 

1.466  1.468 

1.080 

1.479 

1.480 

125.6 

123.0  123.7 

117.3 

109.6 

111.9 

-619.127295 

113.3 

117.2 

108.2 

106.9 

C  Si  Si 

1.702  2.339 

1.076  1.077 

1.473 

1.475 

1.477 

123.0 

122.6  122.7 

119.5 

108.9 

110.5 

-619.121299 

114.7 

117.5 

109.4 

108.3 

Si  Si  Si 

2.133  2.336 

1.469  1.469 

1.472 

1.476 

1.476 

122.7 

122.1  123.8 

118.3 

108.2 

110.8 

-870.165074 

114.2 

119.1 

109.3 

108.3 

Appendix  6.  Geometries  and  Energies  of  Allylic  Radicals 


H 

1 

h2 

Hj\  / 
1  s 

H,  H, 

X  Y  Z 

XY  YZ 

H1X  H2X 

HY 

H1Z  H2Z 

XYZ 

H1XYH2XY 

HYX 

H1ZY  H2ZY 

E(RHF/6-31G*0 

HXH 

HYZ 

HZH 

C  C  C 

1.391  1.391 

1.076  1.074 

1.078 

1.076  1.074 

124.5 

121.2  121.5 

117.7 

121.2  121.5 

-116.468100 

117.4 

117.7 

117.4 

C  Si  C 

1.778  1.778 

1.077  1.076 

1.466 

1.077  1.076 

122.7 

122.2  122.5 

118.7 

122.2  122.5 

-367.477722 

115.2 

118.7 

115.2 

Si  C  C 

1.849  1.353 

1.477  1.476 

1.080 

1.077  1.077 

124.6 

112.2  112.6 

117.8 

122.0  122.0 

-367.496986 

111.4 

117.4 

116.0 

Si  Si  C 

2.281  1.796 

1.477  1.477 

1.475 

1.077  1.077 

119.6 

113.6  113.8 

116.5 

122.2  122.3 

-618.523178 

110.2 

114.0 

114.8 

Si  C  Si 

1.813  1.813 

1.473  1.473 

1.082 

1.473  1.473 

127.2 

114.6  115.3 

115.4 

114.6  115.3 

-618.519739 

112.5 

115.4 

112.5 

Si  Si  Si 

2.287  2.287 

1.477  1.477 

1.479 

1.477  1.477 

119.1 

114.5  114.2 

112.0 

114.5  114.2 

-869.571405 

110.1 

112.0 

110.1 

Appendix  7. 


x  y  z 

E(RHF/6-31G*) 

c  c  c 

-116.438510 

Si  C  C 
-367.440769 

C  Si  C 
-367.455668 

Si  Si  C 
-618.493100 

C  C  Si 
-367.491447 

Si  C  Si 
-618.506179 

C  Si  Si 
-618.506562 

Si  Si  Si 


Geometries  and 


Energies  of  Twisted  Aliylic  Systems 


h2 


H 

I 


H 


H, 


XY  YZ 

H1X  H2X 

HY 

HZ 

XYZ 

H1XYH2XY 

HYX 

HZY 

HXH 

HYZ 

HZH 

1.326  1.479 

1.076  1.075 

1.081 

1.076 

124.8 

121.8  121.4 

115.5 

120.5 

116.8 

119.7 

117.4 

1.700  1.477 

1.467 

1.467 

1.083 

1.076 

127.4 

123.7 

121.4 

116.4 

120.6 

114.9 

116.2 

117.0 

1.693  1.851 

1.076 

1.076 

1.469 

1.077 

125.2 

123.3 

121.6 

120.4 

122.6 

115.1 

114.4 

114.9 

2.126  1.854 

1.468 

1.468 

1.470 

1.076 

125.1 

123.2 

121.3 

119.3 

122.3 

115.5 

115.6 

115.3 

1.325  1.878 

1.076 

1.077 

1.081 

1.479 

124.1 

122.3 

122.1 

117.7 

111.7 

115.6 

118.2 

108.7 

1.695  1.858 

1.466 

1.468 

1.081 

1.480 

125.7 

122.9 

122.7 

117.2 

112.8 

114.4 

117.1 

108.0 

1.701  2.329 

1.076 

1.076 

1.473 

1.477 

122.2 

122.7 

122.5 

119.3 

113.2 

114.8 

118.6 

109.3 

1.468  1.468 

122.2  123.6 

114.2 


-869.550101 


2.132  2.324 

121.7 


1.473 

118.1 

120.2 


1.476 

113.5 

109.4 


Appendix  8.  Geometry  information  for  H,PBHj. 
All  lengths  are  in  A  and  angles  are  in  degrees. 


Geometry 

name 

PB 

PH 

PH’ 

BH 

BH’ 

HPB 

H’PB 

PBH 

PBH’ 

HPBH 

H’PBH 

Ground  State 

1.905 

1.399 

1.187 

102.2 

120.0 

38.8 

RHF  Optimized 

1.399 

1.187 

102.2 

120.0 

141.2 

Ground  State 

1.863 

1.406 

1.191 

104.7 

119.7 

36.8 

MP2  Optimized 

1.406 

1.191 

104.7 

119.7 

143.2 

Planar  Phosph 

1.808 

1.380 

1.184 

124.9 

118.4 

0.0 

RHF  Optimized 

1.380 

1.184 124.9 

118.4 

180.0 

Planar  Phosph 

1.785 

1.393 

1.189 

125.4 

118.1 

0.0 

MP2  Optimized 

1.393 

1.189 

125.4 

118.1 

180.0 

Rotated  Boron 

1.973 

1.409 

1.188 

95.2 

120.3 

47.2 

RHF  Optimize 

1.409 

1.187 

95.2 

120.8 

-47.2 

Rotated  Planar 

1.961 

1.378 

1.185 

122.5 

120.3 

90.0 

RHF  Optimized 

1.378 

1.185 

122.5 

120.3 

-90.0 

Appendix  9.  Geometry  information  for  HP(BH,),. 
All  lengths  are  in  A  and  angles  are  in  degrees. 


Geometry  Name 

PH 

PB 

BH 

HPB 

BPB 

PBH 

HPBH 

BH’ 

PBH’ 

HPBH’ 

Ground  State 

1.391 

1.869 

1.186 

111.9 

116.2 

120.4 

-19.7 

RHF  Optimized 

1.186 

118.3 

164.5 

Ground  State 

1.401 

1.837 

1.190 

114.8 

118.6 

120.8 

-14.4 

MP2  Optimized 

1.190 

117.2 

169.0 

Planar  Phosp 

1.385 

1.849 

1.185 

117.3 

125.5 

120.2 

0.0 

RHF  Optimized 

1.185 

118.0 

180.0 

Planar  Phosp 

1.399 

1.829 

1.189 

117.7 

124.6 

120.6 

0.0 

MP2  Optimized 

1.190 

117.2 

180.0 

Monorotated 

1.404 

1.900 

1.189 

102.1 

105.5 

120.4 

-35.9 

1.187 

119.8 

151.1 

(The  rotated  Boron) 

1.969 

1.187 

119.5 

115.5 

1.187 

121.8 

-63.8 

Planar  Monorot 

1.387 

1.811 

1.185 

120.2 

128.1 

118.7 

0.0 

1.186 

118.8 

180.0 

(The  rotated  Boron) 

1.955 

1.185 

111.7 

120.2 

90.0 

1.185 

120.2 

-90.0 

Dirotated 

1.419 

1.969 

1.188 

91.4 

103.1 

120.0 

124.6 

1.188 

121.1 

-54.9 

Planar  Dirot 

1.384 

1.961 

1.186 

116.2 

127.6 

120.6 

90.0 

1.186 

120.6 

-90.0 

Appendix  10.  Geometry  information  for  PCBH^, 

All  lengths  are  in  A  and  angles  are  in  degrees.  (Note,  B’  is  the  rotated  boron.) 


Geometry  Name  PB 

Ground  state  1.873 
Mono-rot  planar  1.853 


PB’ 

BH 

BPB 

B’H 

1.873 

1.186 

120.0 

1.940 

1.186 

120.5 

1.186 

1.951 

1.185 

117.7 

1.186 

BPB’ 

PBH 

BPBH 

PB’H 

BPB’H 

120.0 

119.4 

0.0 

119.7 

117.9 

0.0 

120.3 

90.0 

117.9 

120.9 

11.3 

120.2 

104.6 

Mono- rot  pyramid  1.859 
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